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1 Introduction 



The "Workshop on Precision Measurements of a^" was held at the Max Planck Institute for 
Physics, February 9th through 11th of 2011. The meeting brought together experts from 
several different fields where high precision measurements of as{mz) are currently being 
made. Main goals of the workshop were to facilitate discussion between the groups, and in 
particular to give speakers the opportunity to explain details that one would normally not 
be able to present at a conference, but which have an important impact on the analyses. In 
each field the session was led off by a review speaker, followed by more specialized talks, and 
was closed with a dedicated time period for discussions. 

There were 67 physicists who took part in the workshop, and 35 talks were presented. 
Slides as well as background reference materials are available on the conference website 

http:/ /www.mpp. mpg.de/alphas 

The sessions and talks in the workshop program were 

• Welcome 

— "World Summary of as (2009) and beyond", S. Bethke 

• Os from Deep Inelastic Scattering and Global Fits 

— "Review of Determinations from Jets at HERA" by C. Glasman 

— "q!s from Deep-Inelastic Scattering: DESY Analysis" by J. Bliimlein 

— "CTEQ-TEA Parton Distribution Functions and a<j" by CP. Yuan 

— "a;^ in MSTW Analyses" , A. Martin 

— "Unbiased from Global Fits: The NNPDF Approach" by S. Forte 

— "Hera PDF" by B. Reisert 

— Discussion Session on DIS and Global Fits, convened by V. Radescu 

• Measurements of Og from r Decays 

— "a^ Determinations from Hadronic r Decays" by A. Pich 

— "a^ from Contour Improved Perturbation Theory (CIPT)" by S. Descotes-Genon 

— "Fixed Order Perturbation Theory (FOPT) Analysis" by M. Beneke 

— "FOPT and CIPT in r Decays" by S. Menke 

— "Duality Violations in Hadronic r Decays" by M. Goltermaan 

— "Perturbative Input to r Decays" by J. Kiihn 

— "Running and Decoupling of as at Low Scales" by M. Steinhauser 

— Discussion Session on r Decays, convened by A. Hocker 

• as from Z Decays and Electroweak Observables 

— "a^ in Electroweak Physics" by J. Kiihn 

— "a^ with Global Analysis of Particle Properties (GAPP)" by J. Erler 

— "a^ from the Hadronic Width of the Z" by K. Monig 

— Discussion Session on Elecroweak Analyses, convened by W. HoUik 

• as from Event Shape Measurements 

— "Review of event-shape measurements of Os" by G. Salam 

— "a^ at NNLO and NNLA from (mainly) ALEPH data" by T. Gehrmann 

— "NNLO and Classic Power Corrections" by B. Webber 

— "q!s from Soft-Collinear Effective Theory analysis" by V. Mateu 
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— "Five jets at LEP at NLO and asimz)" by G. Zanderighi 

— "Jet and Event Shape Observables at the LHC" by K. Rabbertz 

— "Monte Carlo Tuning with Professor" by H. Schulz 

— "Experimental Issues and Combination of Results" by S. Kluth 

— Discussion Session on Event Shapes, convened by Y. Dokshitzer 

— MPI Colloquium, "On the Other Side of Asymptotic Freedom" by Y. Dokshitzer 

• Lattice Calculations and Quarkonia 

— "Lattice QCD Calculations and a^" by A. Kronfeld 

— "The QCD Coupling from Lattice QCD (HPQCD)" by CP. Lepage 

— "a, from Lattice QCD (PACS-CS)" by S. Aoki 

— "a, from Lattice QCD (JLQCD)" by S. Hashimoto 

— "a, from Lattice QCD (ALPHA)" by R. Sommer 

— "a^ from Bottomonium" by N. Brambilla 

— Discussion Session on Lattice QCD, convened by C. Davies 

• Applications and Future of Measurements 

— "Os for New Physics" by Y. Nomura 

— "The Strong Coupling and LHC Cross Sections" by F. Petriello 

— "cts at the International Linear Collider (ILC)" by C. Pahl 

— "The Principle of Maximum Conformality" by S. Brodsky 

— Final Discussion Session, convened by S. Bethke 

This web proceedings represent a collection of extended abstracts and references for the 
presentations, summarizing the most important results and issues. In these writeups, unless 
otherwise indicated, the strong coupling is defined in the MS scheme, and as{mz) is 

the coupling with five light quark flavors. 
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Workshop on Precision Measurements as 2011: 
Welcome and Introduction 

Siegfried Bethke 
Max-Planck-Institute of Physics, Munich, Germany 

Heartiest welcome to all participants of the 2011 Workshop on Precision Measurements 
on as at the Max-Planck-Institute of Physics at Munich! The principal organisers, i.e. Andre 
Hoang (Vienna Univ.), Stefan Kluth (MPP Munich), Jochen Schieck (LMU Munich), Iain 
Stewart (MIT) and myself, are very glad to see so many representatives and world-experts 
from all subfields of the workshop topics present at this meeting. 

Measurements of a^, one of the basic constants of nature whose numerical value is not 
given by any theory, continue to be at the fore-front of theoretical and experimental endeavors 
in particle physics. The purpose of this workshop is 

• to review the latest theoretical developments and experimental studies in this field, 

• to discuss and possibly solve open questions and issues which currently limit the pre- 
cision to which Us is known, 

• to provide input for calculating a new 2011 world average of as{mz) 

The latest comprehensive summary of as measurements and calculation of a world com- 
bined average value of asijnz) was given in P in 2009, which resulted in 

as{mz) = 0.1184 ±0.0007 . 

This result was also adopted by the latest review of the particle data group [2]. 

Table 1: Summary of recent measurements of as{mz)- The rightmost two columns give the 
exclusive mean value of as{mz) calculated without that particular measurement, and the number 
of standard deviations between this measurement and the respective exclusive mean. 



Process 


Q [GeV] 


Osimz) 


excl. mean as{mz) 


std. dev. 


r-decays 


1.78 


0.1197 ±0.0016 


0.11818 ±0.00070 


0.9 


DIS [F2] 


2 - 170 


0.1142 ±0.0023 


0.11876 ±0.00123 


1.7 


DIS [e-p — )• jets] 


6 - 100 


0.1198 ±0.0032 


0.11836 ±0.00069 


0.4 


Lattice QCD 


7.5 


0.1183 ±0.0008 


0.11862 ±0.00114 


0.2 


T decays 


9.46 




0.11841 ± 0.00070 


0.1 


e"'"e~ [jets & shps] 


14 - 44 


0.1172 ±0.0051 


0.11844 ± 0.00076 


0.2 


e'''e~ [ew prec. data] 


91.2 


0.1193 ±0.0028 


0.11837 ±0.00076 


0.3 


e"'"e~ [jets & shps] 


91 - 208 


0.1224 ±0.0039 


0.11831 ±0.00091 


1.0 
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Figure 1: Summary of measurements of as{mz)- The vertical line and shaded band mark the 
final world average value of as{mz) = 0.1184 ±0.0007 determined from these measurements. 
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The individual results entering that determination of as{mz) are given in Table 1, and a 
graphical representation in Fig. 1; see [T] and references therein. All results are compatible 
with the quoted overall world average; the largest deviation from an exclusive mean being 
less than two standard deviations (interpreting quoted uncertainties of results as gaussian 
errors). In fact, the of the overall average is close to unity per degree of freedom [1]. While 
the 2009 average came out in a rather consistent way, some questions and potential problems 
remained: 

• the individual result with smallest quoted uncertainties is from lattice QCD; it largely 
dominates (not so much) the overall average and (mainly) its assigned error; 

• there are large systematic differences between different studies of as from r decays, 
depending on the type of perturbative prediction (CIPT and FOPT), leading to an 
increased uncertainty assigned to the value listed in Table 1; 

• are there systematic unknowns in the DIS/F2 result? 

• theoretical uncertainties are defined using different definitions in most of the input 
studies; do they lead to a consistent guess-timate of the overall uncertainties? 

As will be seen in the course of this workshop, the consistency of an overall new world 
summary may look more unfavorable in the light of some of the most recent results (see the 
summary contribution to this workshop [3]). 

References 

[1] S. Bethke, Eur.Phys.J. C64 (2009) 689-703, arXiv:0908.1135 [hep-ph]. 
[2] Particle Data Group, J.Phys.G G37 (2010) 075021. 
[3] S. Bethke, these proceedings. 
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Review of ag determinations from jets at HERA 



Claudia Glasman 
(On behalf of the HI and ZEUS Collaborations) 

Universidad Autonoma de Madrid, Spain 

The HI and ZEUS Collaborations at HERA have performed precise determinations of 
as{Mz) using observables such as jet cross sections, ratios of jet cross sections, and internal 
structure of jets in different regimes. Some of them are reviewed in this report. In addition, 
the energy-scale dependence of the coupling has been determined very precisely. The results 
are in good agreement with the predicted running of with small experimental uncertainties 
in a wide range of the scale. 

Values of as{Mz) were extracted by the HI Collaboration from the double-differential 
inclusive-jet, dijet and trijet cross sections at low (5 < < 100 GeV^) using 43.5 pb~^ 
of integrated luminosity [1] . Jets were searched using the kx cluster algorithm in the Breit 
frame with Pt > 5 GeV and —1 < ?7lab < 2.5. The values of as{Mz) extracted are shown 
in Table 1 (rows 1, 2 and 3; row 4 shows the combined value). These values have small 
experimental uncertainties; however, the theoretical uncertainty, dominated by terms beyond 
NLO, is large. A reduction of the theoretical uncertainties can be achieved by determining 
as{Mz) from the trijet to dijet ratio; the value obtained is shown in row 5 and has a smaller 
theoretical uncertainty. 

Values of as{Mz) were extracted by the HI Collaboration from the double-differential 
inclusive-jet, dijet and trijet normalised cross sections at medium (150 < < 15000 GeV^) 
using 395 pb~^ of integrated luminosity [2]. Jets were searched using the cluster algo- 
rithm in the Breit frame with Pt > 5 GeV and —0.8 < ?7lab ^ 2. The values of as{Mz) 
extracted are shown in rows 6, 7 and 8. Row 9 shows the combined as{Mz) value, which has 
a very small experimental uncertainty; the theoretical uncertainties, though still dominated 
by the terms beyond NLO, are smaller than for the low analysis. This reduction was 
accomplished by using normalised cross sections, for which correlated uncertainties cancel, 
and making the extraction at higher Q^, where the contribution from the terms beyond NLO 
is reduced. 

The ZEUS Collaboration extracted values of as{Mz) from the inclusive-jet cross sections 
based on the /c^, anti-Zc^ and SIScone jet algorithms at high (Q^ > 500 GeV^) using 
82 pb~^ of integrated luminosity [3],|1]. Jets were searched in the Breit frame with -E^g > 8 
GeV and —2 < rf^^ < 1.5. The values of as{Mz) extracted are shown in rows 10, 11 and 
12. Row 13 shows an updated value using HERA II statistics {C = 300 pb-^) [5]. The 
experimental uncertainty is dominated by the jet energy scale and amounts to 2%. The 
theoretical uncertainties are dominated by the terms beyond NLO and amount to around 
1.5% for the three jet algorithms. The reduction of the theoretical uncertainty was obtained 
by restricting to the high region; even though the experimental uncertainty increases at 
high (mainly due to statistics), the total uncertainty remains small. 

Values of as{Mz) were extracted by the ZEUS Collaboration from the inclusive-jet cross 
section in photoproduction at high E^r^^ (21 < i?;^* < 71 GeV) using 189 pb~^ of integrated 
luminosity [6], [7]. Jets were searched using the krp, anti-fcr and SIScone jet algorithms in the 
laboratory frame. The values of as{Mz) extracted are shown in rows 14, 15 and 16. These 
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have similar precision. In photoproduction, there is an additional uncertainty coming from 
the photon PDFs, which is as large as that coming from higher orders. 

A completely different approach to the extraction of at HERA is given by the analysis 
of the internal structure of jets. The internal structure of jets can be studied by means of the 
integrated jet shape and the subjet multiplicity. The ZEUS Collaboration extracted values of 
as{Mz) from the measured mean integrated jet shape [S] for E-!^^ > 21 GeV using 82.2 pb~^ 
of integrated luminosity and from the mean subjet multiplicity [9] for i?;^* > 25 GeV with 
38.6 pb~^ of integrated luminosity. Jets were searched using the kx cluster algorithm in 
the laboratory frame. The values of as{Mz) extracted are shown in rows 17 and 18. In 
both cases, the experimental uncertainties are small; however, the theoretical uncertainty, 
dominated by the terms beyond NLO, is much bigger than for the extraction from jet cross 
sections. 

From the results presented, it is concluded that the dominant uncertainty in the extrac- 
tion of as{Mz) at HERA comes from terms beyond NLO. This uncertainty decreases with 
increasing or so an extraction at high or i?^* is advantageous to minimise this 
contribution. In addition, such uncertainty cancels partially when ratios of jet cross sections 
are used for the extraction, and lower or Ei^^ values can be used. The uncertainties 
coming from the PDFs are also smaller as increases. 

To summarise, jet physics at HERA provides precise values of as{Mz) in different regimes 
and precise determinations of the running of the coupling over a wide range of the scale. 
There is still room for improvement. All HERA data have been analysed: the values of 
as{Mz) obtained have very small experimental uncertainties; the theoretical uncertainties 
are dominant, in particular that coming from terms beyond NLO. Therefore, enormous 
benefit will be gained from NNLO calculations for jet cross sections at HERA to improve 
the precision in as{Mz). 

References 

[1] HI Collaboration, F.D. Aaron et al., Eur. Phys. J. C 67 (2010) 1. 

[2] HI Collaboration, F.D. Aaron et al., Eur. Phys. J. C 65 (2010) 363. 

[3] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 649 (2007) 12. 

[4] ZEUS Collaboration, H. Abramowicz et al., Phys. Lett. B 691 (2010) 127. 

[5] ZEUS Collaboration, H. Abramowicz et al., 

http: / /www-zeus. desy.de/public_results/functiondb. php?id=ZEUS-prel- 10-002. 

[6] ZEUS Collaboration, H. Abramowicz et al., 

http:/ /www-zeus. desy.de/public_results/functiondb. php?id=ZEUS-prel-10-003. 

[7] ZEUS Collaboration, H. Abramowicz et al., 

http: / / www-zeus. desy. de / public_results/ functiondb .php?id=ZEUS-prel- 10-0 15 . 

[8] ZEUS Collaboration, S. Chekanov et al., Nucl. Phys. B 700 (2004) 3. 
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[9] ZEUS Collaboration, S. Chekanov et al., Phys. Lett. B 558 (2003) 41. 
[10] C. Glasman, hep-ex/0506035. 

[11] HI and ZEUS Collaborations, Hlprelim-07-132 and ZEUS-prel-07-025, 

http:/ /www-zeus. desy.de/public_results/functiondb. plip?id=ZEUS-prel-07-025. 



Process 


Collab. 


Value 


Exp. 


Th. 


Total 


(%) 


(1) Inc. jets at low 


HI 


0.1180 


0.0018 


+0.0124 
—0.0093 


+0.0125 
—0.0095 


+10.6 
—8.1 


(2) Dijets at low Q"^ 


HI 


0.1155 


0.0018 


+0.0124 
—0.0093 


+0.0125 
—0.0095 


+10.8 
—8.2 


(3) Trijets at low Q'^ 


HI 


0.1170 


0.0017 


+0.0091 
—0.0073 


+0.0093 
—0.0075 


+7.9 
—6.4 


(4) Combined low 


HI 


0.1160 


0.0014 


+0.0094 
—0.0079 


+0.0095 
—0.0080 


+8.2 
—6.9 


(5) Trijet/dijet at low 


HI 


0.1215 


0.0032 


+0.0067 
—0.0059 


+0.0074 
—0.0067 


+6.1 
—5.5 


(6) Inc. jets at medium 


HI 


0.1195 


0.0010 


+0.0052 
—0.0040 


+0.0053 
—0.0041 


+4.4 
—3.4 


(7) Dijets at medium 


HI 


0.1155 


0.0009 


+0.0045 
—0.0035 


+0.0046 
—0.0036 


+4.0 
—3.1 


(8) Trijets at medium 


HI 


0.1172 


0.0013 


+0.0053 
—0.0032 


+0.0055 
—0.0035 


+4.7 
—3.0 


(9) Combined medium 


HI 


0.1168 


0.0007 


+0.0049 
—0.0034 


+0.0049 
—0.0035 


+4.2 
—3.0 


(10) Inc. jets at high (anti-A;^) 


ZEUS 


0.1188 


+0.0036 
—0.0035 


+0.0022 

— U.UUZZ 


+0.0042 
—0.0041 


+3.5 
—3.5 


(11) Inc. jets at high Q"^ (SIScone) 


ZEUS 


0.1186 


+0.0036 
-0.0035 


+0.0025 
-0.0025 


+0.0044 
-0.0043 


+3.7 
-3.6 


(12) Inc. jets at high (^^. HERA I) 


ZEUS 


0.1207 


+0.0038 
-0.0036 


+0.0022 
-0.0023 


+0.0044 
-0.0043 


+3.6 
-3.6 


(13) Inc. jets at high (^^. HERA II) 


ZEUS 


0.1208 


+0.0037 
-0.0032 


+0.0022 
-0.0022 


+0.0043 
-0.0039 


+3.6 
-3.2 


(14) Inc. jets in 7p (anti-A;T) 


ZEUS 


0.1200 


+0.0024 
-0.0023 


+0.0043 
-0.0032 


+0.0049 
-0.0039 


+4.1 
-3.3 


(15) Inc. jets in 7p (SIScone) 


ZEUS 


0.1199 


+0.0022 
-0.0022 


+0.0047 
-0.0042 


+0.0052 
-0.0047 


+4.3 
-3.9 


(16) Inc. jets in 7p (/ct) 


ZEUS 


0.1208 


+0.0024 
-0.0023 


+0.0044 
-0.0033 


+0.0050 
-0.0040 


+4.1 
-3.3 


(17) Jet shape 


ZEUS 


0.1176 


+0.0013 
-0.0028 


+0.0091 
-0.0072 


+0.0092 
-0.0077 


+7.8 
-6.5 


(18) Subjet multiphcity 


ZEUS 


0.1187 


+0.0029 
-0.0019 


+0.0093 
-0.0076 


+0.0097 
-0.0078 


+8.2 
-6.6 


HERA average 2004 




0.1186 


±0.0011 


±0.0050 


±0.0051 


±4.3 


HERA average 2007 




0.1198 


±0.0019 


±0.0026 


±0.0032 


±2.7 



Table 1: Values of as{Mz) extracted from jet observables at HERA together with their 
uncertainties (rows 1 to 18). The 2004 [10] and 2007 [TT] HERA averages are shown in the 
last two rows. 
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q:s(M|) in NNLO Analyses of Deep-Inelastic World Data 

S. Alekhin^'2, J. Bliimlein \ H. B6ttcher\ and S.-O. Moch^ 

^ Deutsches Elektronen-Synchrotron, DESY, Platanenallee 6, D- 15738 Zeuthen, Germany 
^ Institute for High Energy Physics, 142281 Protvino, Moscow Region, Russia 

The present world data of deep-inelastic scattering (DIS) reached a precision which allows the 
measurement of as(M|) from their scaling violations with an error of (5as(M|) ^ 1%. This 
requires at least NNLO analyses, since NLO fits exhibit scale uncertainties of Arjas(M|) ~ 
0.0050. The NNLO values for obtained are summarized in the following Table. 









BBG 


0.1134 


valence analysis, NNLO [T] 


GRS 


0.112 


valence analysis, NNLO [2] 


ABKM 


0.1135 ±0.0014 


HQ: FFNS A^/ = 3 [3] 


ABKM 


0.1129 ±0.0014 


HQ: BSMN-approach (3] 


JR 


0.1124 ±0.0020 


dynamical approach (3] 


JR 


0.1158 ±0.0035 


standard fit [1] 


MSTW 


0.1171 ±0.0014 


m 


ABM 


0.1147 ±0.0012 


FFNS, incl. combined Hl/ZEUS data [6] 


Gehrmann et al. 


0.1153 ±0.0017 ±0.0023 


e+e~ thrust [7] 


Abbate et al. 


0.1135 ±0.0011 ±0.0006 


e^e~ thrust [8] 


BBG 


n 1141 +00020 
"•-■^^^^ -0.0022 


valence analysis, N'^LO [1] 


world average 


0.1184 ±0.0007 


m 



NNLO non-singlet data analyses have been performed in [T] and [2]. The analysis of Ref. [T] 
is based on an experimental combination of flavor non-singlet data referring to F^^'^lx, Q"^) 
for X < 0.35 and using the respective valence approximations for x > 0.35. The d — u 
distributions and the 0{a1) heavy flavor corrections were accounted for. At low and at 
large x also at low higher twist corrections have to be taken into account [lO]. The 
corresponding region was cut out in [1] performing the fits for the leading twist terms only. 
The analysis could be extended to N^LO effectively due to the dominance of the Wilson 
coefficient in this order [TT] if compared to the anomalous dimension, cf. [1] and [12]. This 
analysis led to an increase of a^(M|) by ±0.0007 if compared to the NNLO value. 

A combined singlet and non-singlet NNLO analysis based on the DIS world data, in- 
cluding the Drell-Yan and di-muon data, needed for a correct description of the sea-quark 
densities, was performed in [3]. In the fixed fiavor number scheme (FFNS) the value of 
as(M|) is the same as in the non-singlet case [1]. The comparison between the FFNS and 
the BMSN scheme [13] for the description of the heavy fiavor contributions induces a sys- 
tematic uncertainty Aa<i(M|) = 0.0006. The NNLO analyses of Ref. [1] are statistically 
compatible with the results of [1] and [3] while those of [5] yield a higher value. 

In Ref. [6] the combined HI and ZEUS data were accounted for in a NNLO analysis for 
the first time, which led to a shift of ±0.0012. However, running quark mass effects [T3| and 
the account of recent data reduce this value again to the NNLO value given in [3]. We 
mention that other recent NNLO analyses of precision data, as the measurement of as(M|) 
using thrust in high energy e^e~ annihilation data [7], [8] result in lower values than the 
2009 world average [9] based on NLO, NNLO and N^LO results. The sensitivity of the fits 
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to a precise description of the longitudinal structure function Fl has been demonstrated in 
[T5] recently, in the case of the NMC data. Inconsistent descriptions of induce a high 
value of as of ~ 0.1170 to be compared with that obtained in [5]. It is observed that the 
values of 



a,(M|) 


with CTNMC 


with F^^^v^^ 


difference 


NNLO 
NNLO +Fl 0(a3) 


0.1135(14) 
0.1122(14) 


0.1170(15) 
0.1171(14) 


+0.0035 ~ 2.3o- 
+0.0050 ~ 3.6a 



as found in NLO fits are systematically higher than those in NNLO analyses, ag mea- 
surements based on jet data can be performed presently at NLO only. Here typical values 
obtained are as(M|) = 0.1156 taoow [IS]; '^s(^l) = 0.1161 ^l^o^ooS [IZ] in recent examples. 
The precise knowledge of as(M|) is of instrumental importance for the correct prediction of 
the Higgs boson cross section at Tevatron and the LHC |18j . 
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Abstract 

After summarizing a simple method to fully reproduce the correlated dependence of 
theoretical cross sections on the QCD coupling strength and parameters of parton 
distribution functions (PDFs), we present a series of CTEQ6.6AS and CTIO.AS PDFs, 
realizing this approach, for a<j values in the interval 0.113 < as{Mz) < 0.123. If is 
solely determined by the QCD global analysis data, and not by the world average value 
of as{Mz) as done in the CTIO PDF analysis, we find that a{Mz) = 0.1197 ± 0.0061 
at the 90% CL. 

In a recent paper [1], we examined the dependence of parton distribution functions (PDFs) 
on the value of the QCD coupling strength as{Mz), in the QCD global analysis performed by 
the CTEQ-TEA (Tung et al.) group. We demonstrated a simple method that is rigorously 
valid in the quadratic approximation normally applied in PDF fitting, and fully reproduces 
the correlated dependence of theoretical cross sections on and PDF parameters. This 
method is based on a statistical relation that allows one to add the uncertainty produced by 
as, computed with some special PDF sets, in quadrature with the PDF uncertainty obtained 
for a fixed as value (such as the CTEQ6.6 PDF set). A series of four CTEQ6.6AS PDFs 
realizing this approach, for as values between 0.116 and 0.120, was made available online 
[2] . Using these PDFs, the combined as and PDF uncertainty can be assessed for theoretical 
predictions at the Fermilab Tevatron and Large Hadron Collider. In a more recent paper 
[3], we provided a similar series, CTIO.AS PDFs, but with a larger range of as values, from 
0.113 to 0.123 [g. 

Below, we summarize a few important features of, for example, the CTEQ6.6AS PDF 
analysis. Using the setup of the CTEQ6.6 PDF analysis and taking the world-average (w.a.) 
value [5] (tts)^^ ± (^tts)^a = 0.118 ± 0.002 (at 90% CL [6].) as a separate data point in 
addition to the complete CTEQ6.6 set of hadronic scattering data, we allow as{Mz) to vary 
within the global fits and find 

asiMz) = 0.1180 ± 0.0019 (at 90% CL). (1) 

Thus, the constraint on as{Mz) is dominated by the world-average uncertainty, {Sas)w.a.- 
Using the ag series of PDFs, and including the PDF uncertainty, we can estimate the uncer- 
tainties of cross section calculations. For any calculated quantity cr, we denote the central 
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prediction, corresponding to as{Mz) = 0.118, by (Tq. There are two contributions to the un- 
certainty: the PDF uncertainty (Actpdf ) and the as uncertainty {AaaJ of a. The combined 
uncertainty Aa for CTEQ6.6+CTEQ6.6AS is {Aaf = (Actpdf)^ + (Aa^f ■ 

The independence of the uncertainty from the PDF uncertainty in the CTEQ6.6AS 
method does not preclude existence of some correlation between the and PDF parameters. 
This correlation arises from the hadronic scattering experiments, which probe a variety of 
combinations of the PDFs and a^, and can be closely examined by analyzing the correlation 
cosine as introduced in Refs. [7]. For that purpose, we performed the full fit CTEQ6.6FAS 
with the floating as and examined the correlation between as{Mz) and individual PDFs 
fa{x,Q). Fig. 8 of [I] shows the correlation cosine, cosy? versus x, for the PDFs that 
have the largest correlations with as{Mz), at Q = 2 and 85 GeV. The best-fit value of 
as{Mz) is thus determined by several types of the data, probing the gluon evolution in 
DIS at moderately small x, the singlet PDF evolution in DIS at large x, and HERA charm 
semi-inclusive DIS data. The correlation of each kind disappears if the relevant data set is 
removed. 

After the completion of Refs. [T] and [3], we performed a similar analysis on the deter- 
mination of as{Mz) value but with the CTIO setup. We find that, without including the 
world-average as{Mz) value as an input data, the constraints of the global fit on as are 
relatively weak, and 

as{Mz) = 0.1197 ± 0.0061 (at 90% CL), (2) 

where the error is an estimate of the uncertainty (at about 90 % CL.) from all typical 
sources, including the statistical error, dependence on the parametrization form and free 
theoretical parameters. Within this range, x2 exhibits nearly parabolic dependence on as- 
Eq. [2] shows clearly that constraints imposed on as by the hadronic scattering data only are 
significantly weaker than from the world average data point. 

Finally, we note that replacing the NMC F2 data by the corresponding differential cross 
section data in the CTIO. AS PDF global analysis does not change our conclusion about the 
determination of the as{Mz) value. 
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In the MSTW2008 P global PDF analysis of DIS and related hard-scattering data the 
value of as{M'^) was left as a free parameter and its optimum value determined. Subse- 
quently, in [2], a detailed study was made of its 'experimental' uncertainty. The values 
found in the NLO and NNLO analyses are 

NLO: as{Ml) = 0.1202 (68% C.L.) (90% C.L.), (1) 

NNLO: a5(M|) = 0.1171 +Hoi4 (68% C.L.) ^^osl (90% C.L.). (2) 

We note that the NNLO value is smaller than that at NLO since NNLO evolution is faster 
and since NNLO coefficient functions are generally positive. Fig. [1] shows the constraints 
on as coming from the individual data sets in the NNLO global analysis. Inclusive jet 
production is the only process proportional to a| at LO. Moreover, since the normalisation 
of the CDF jet production data is tied to the observed CDF Z rapidity distribution, we 
see that these jet data give the tightest constraint on as- In addition to the 'experimental' 
errors shown in ([T]) and ([2]), there is also a 'theory' uncertainty, which is estimated [2] to be 
±0.003 at NLO and at most ±0.002 at NNLO. 

The as{M^) values found by MSTW are, in general, greater than those of other PDF 
analyses, and are more consistent with the world average value [3]. There are two reasons 
for this. First, MSTW have a more flexible low-x parametrisation of the input gluon PDF. 
The flexibility is required by the data and gives a negative input behaviour at small x, which 
is conflrmed by the NNPDF analyses The second reason is the inclusion of the Tevatron 
jet data in the analysis. It is informative to repeat the NNLO MSTW2008 analysis fltting 
only inclusive DIS data, which gives asiM"^) = 0.1104 compared to asiM"^) = 0.1171 for 
the global flt, but with a negative input gluon for x > 0.4 due to lack of a data constraint, 
implying a negative F^^'^'^'^ and an awful description of Tevatron jet data (x^/^pts. ~ 10). 
Fixing the high-x gluon parameters to be the same as in the global flt gives a5(M|) = 0.1172 
in the DIS-only flt. A DIS-only flt without BCDMS data gives as(M|) = 0.1193, while a 
global flt without BCDMS data gives as{M^) = 0.1181. We conclude that the Tevatron jet 
data are vital to pin down the high-x gluon, giving a smaller low-x gluon by the momentum 
sum rule, at the expense of some deterioration in the flt quality of BCDMS data. Thus, both 
reasons imply a larger as from the scaling violations (~ asg) of HERA data. 

Concerning the NMC data, ref. [5] notes that it is better to use the NMC cross-section 
measurements rather than F2, since for x < 0.12 NMC used a Q^-independent R = (Tl/ut 
value to extract their default F2. Does this bias the MSTW value of a^? As a check, we 
repeated the NNLO global flt with the NMC F2 extracted using the SLAG i?i99o(x, Q^) for 
all X, which is close to the MSTW NNLO i?(x, Q"^) in the most relevant x bins. This has a 
very small effect, with the NNLO a5(M|) moving only from 0.1171 to 0.1167. 

Finally, let us comment on the 'low' value as{Ml) = 0.1141+|];[^[522 

obtained in [6], since 

it is often taken as the 'DIS' value in world average compilations, where it is listed as DIS 
F2(N^L0) [5]. This determination assumes that F2 is pure non-singlet for x > 0.3. We may 
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Figure 1: Ranges of as(M|) for which data sets are described within their 90% C.L. limit (outer 
error bars) or 68% C.L. limit (inner error bars) in the NNLO global fit. The points (•) indicate 
the values of asiM^) favoured by each individual data set n, that is, the values for which Xn 
minimised. The figure is taken from [2]. 

examine this assumption using MSTW2008 PDFs. For the F| and measurements with 
X > 0.3 we have = 329 for 282 data points, of which 160 are from BCDMS. However, 
if we consider only non-singlet contributions then = 1449 for 282 data points. In fact, 
contributions other than valence quarks amount to about 10% at a; = 0.3, and still 2% at 
X = 0.5. We conclude that the anomalously low value of as found in [6] is due both to the 
dominance of BCDMS data (cf. Fig. [1]) and to the neglect of singlet contributions. 
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The determination of as (Mz) from the same wide of set of data which is used to determine 
PDFs is appealing because it simultaneously exploits the dependence on the coupling of 
scaling violations as well as that on individual hard matrix elements of the various processes 
under consideration. On the other hand, in such a determination the value of as is necessarily 
correlated to the best-fit form of the PDFs, and thus subject to potential sources of bias 
which may affect the PDFs, such as for example an insufficiently flexible parametrization. 

Determining with NNPDF partons has the advantage that in this approach parametriza- 
tion bias is reduced to a minimum, thanks to the use of neural networks as unbiased interpo- 
lating functions coupled with a Monte Carlo approach. It has the disadvantage that, because 
NNPDF partons are delivered in the form of a Monte Carlo sample, the of the comparison 
between data and theory is a random variable, which only tends to a constant value in the 
limit in which the size of the Monte Carlo sample tends to infinity. In order to determine as 
accurately one thus needs a number of Monte Carlo replicas of the same order of magnitude 
of the number of datapoints (more than 3500 for NNPDF and other global PDF sets). 

We have determined as by constructing sets of at least 500 replicas for each of the 
eleven values of as for which NNPDF2.1 parton distributions have been provided [1]. The 
comparison to the data is done at NLO in QCD and includes heavy quark mass effects using 
the so-called FONLL method [2]. The uncertainty on each value of due to the finite size 
of the replica sample is determined using the bootstrap method from the fiuctuations of 
computed from subsets of replicas in the given sample. The value of as is then determined 
by fitting a parabola to the viewed as a function of as- The minimum of the parabola 
provides the best-fit value of as (Mz) while the Ax^ = 1 range gives the uncertainty on it. 
The further uncertainty due to the finite-size fiuctuations is propagated from the uncertainty 
on each data point and is essentially negligible with the given sample size. The quality of 
the parabolic fit is assessed by evaluating the corresponding XparZ-^dof, with A'dof = —3. 

We have performed [3] this determination both for the global NNPDF2.1 dataset, for all 
deep-inelastic data included in this dataset, and for HERA data only. The best-fit values and 
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as (Mz) 


NNPDF2.1 


0.1191 ± 0.0006*^*^* ± O.OOOP™" 


NNPDF2.1 DIS-only 
NNPDF2.1 HERA-only 


0.1177 ± 0.0009*^*^* ± 0.0002?"^°= 
0.1103 ± O.OOSS'^*^* ± 0.0003?"^°" 



Table 1: Values of (Mz) and associated uncertainties. All uncertainties shown are 68% confi- 
dence levels. 

uncertainties are collected in Table [U the uncertainty from the Ax^ = 1 range denoted with 
"stat" and the uncertainty from the finite size of the replica sample denoted with "proc" . 

Our best-fit value of as is in good agreement with the current PDG value, and it has 
a surprisingly small statistical uncertainty. The statistical uncertainty increases as the size 
of the dataset is reduced, as it ought to. We see no evidence that DIS data prefer a sig- 
nificantly lower value of as'- the difference between the global and DIS-only determinations 
is compatible with a statistical fluctuation. HERA data, however, lead to a rather smaller 
value of as, though with much larger uncertainty. This may be related to deviations from 
predicted scaling violations, which have been observed [4] in HERA data at low x and Q^. 
If so, the value of as could be brought in line with other determinations by an inclusion of 
small X resummation [3]. 

From the behaviour of the as a function of as for individual datasets for the global and 
the DIS-only fit we find evidence [3] that BCDMS data prefer a lower value of as, but only 
in the DIS only fit. Also, the correlation between the of the fit to various datasets and 
individual PDFs as a function of as provides evidence that this behaviour is due to the fact 
that at low the for BCDMS can be lowered by modifying the gluon distribution in a 
way which is incompatible with jet data. However, the best-fit as for the DIS fit is unaffected 
by all this, perhaps becuase of the great flexibility of the NNPDF parton parametrization. 

Theoretical uncertainties on our results are likely to be dominant. The main ones are 
likely to be related to higher order QCD corrections and their resummation, and to the 
treatment of heavy quarks. These uncertainties are presumably of similar size here as in 
other determinations of as based on the same QCD processes. A detailed study of these un- 
certainties, which are likely to be dominant, would be of great interest once QCD corrections 
to the highest available order have been included. 
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The inclusive character of the total r hadronic width renders possible [1] an accurate 
calculation of the ratio Rr = ^[t~ — )■ z/,- hadrons]/r[r~ — i- Ure'Ue]- Its Cabibbo-allowed 
component can be written as [2] 

Rr,v+A = Nc iKdP ^Ew {1 + 5p + 5np} , (1) 

where Nc = 3 is the number of quark colours and Sew = 1.0201 ± 0.0003 contains the 
electroweak radiative corrections. The non-perturbative contributions are suppressed by six 
powers of the r mass [1] and can be extracted from the invariant-mass distribution of the 
final hadrons [3]. From the ALEPH data, one obtains 5np = -0.0059 ± 0.0014 [1]. 
The dominant correction (~ 20%) is the perturbative QCD contribution [1] [3] 

5p = = J2{K^ + gn)a^^ = ^r„<, (2) 

n=l n=l n=l 

which is determined by the coefficients of the perturbative expansion of the {Np = 3) QCD 
Adler function, already known to O(a^) [5]: Kq = Ki = 1; K2 = 1.63982; KsCMS) = 6.37101 
and Ki(MS) = 49.07570. The functions [3] 

-^f -( f 1 - 24 + 4 - 4) = < + o(K») (3) 

2m J\s\=m% s \ n J \ mi mP^ m% J 

are contour integrals in the complex plane, which only depend on a-r = as{frtl)/TT. Using the 
exact solution (up to unknown /3„>4 contributions) for as{—s) given by the renormalization- 
group /3- function equation, they can be numerically computed with very high accuracy [5]. 

If the integrals A'^^\as) are expanded in powers of a^, one recovers the naive perturbative 
expansion of 5p shown in the rhs of Eq. ([2]). This approximation is known as fixed-order 
perturbation theory (FOPT), while the improved expression, keeping the non-expanded values 
of y4*^")(as), is usually called contour-improved perturbation theory (CIPT) [3]. Even at C(a^), 
FOPT gives a rather bad approximation to the integrals ^'^"^(a^), overestimating 5-p by 12% 
at At- = 0.11. The long running of as{—s) along the circle \s\ = generates very large Qn 
coefficients, which depend on Km<n and /3m<n [3]: Qi = 0, g2 = 3.56, = 19.99, (74 = 78.00, 
g5 = 307.78. These corrections are much larger than the original Kn contributions, giving 
rise to a badly behaved perturbative series (at the four-loop level the expansion of as{—s) 
in powers of a,- is only convergent for < 0.11, which is very close to the physical value 
of a^). Thus, it seems compulsory to resum the large logarithms, log" {—s/mi), using the 
renormalization group. This is precisely what CIPT does. 

It has been argued that in the asymptotic regime (large n) the renormalonic behaviour of 
the Kn coefficients could induce cancelations with the running gn corrections, which would 
be missed by CIPT. In that case, FOPT could approach faster the 'true' result provided by 
the Borel summation of the full renormalon series. This happens actually in the large-/3i 
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limit, which however does not approximate well the known perturbative series (for n < 4 the 
true Kn coefficients add constructively with the Qn contributions). Models of higher-order 
corrections which assume a precocious asymptotic behaviour of the Adler function already 
at n = 3,4 [H] [7] seem to favour the FOPT result. The CIPT procedure is much more 
reliable in all other scenarios. 

The present experimental value Rr,v+A = 3.4771±0.0084 [S] implies 6p = 0.2030±0.0033. 
The two different treatments of the perturbative series result in 

a.(m^)ciPT = 0.3412 ± 0.0041,^ 1°:°°^^^^+^^^?^!°:^^^^^^^ = 0.344 ±0.014, (4) 

«.K)poPT = 0.3194 ± 0.0028,^ = 0.321 ±0.015. (5) 

Higher-order corrections have been estimated adding the fifth-order term A^^\as) with 
= 275 ± 400. We have also included the 5-loop variation with changes of the renor- 
malization scale in the range s) G [0.4,2.0]. The error induced by the truncation of 

the /3 function at fourth order has been conservatively estimated through the variation of 
the results at five loops, assuming (3^ = ±f3l/f33 = ±443; in CIPT this slightly changes the 
values of A^^'^as), while in FOPT it increases the scale sensitivity. The FOPT result shows 
as expected [2] [S] a much more sizeable /i dependence, but it gets smaller errors from 6p 
and K^. The three theoretical uncertainties {K^, n, (5^) have been added linearly and their 
sum combined in quadrature with the 'experimental' error from 5p. 

Combining the two results with the PDG prescription (scale factor S = 1.14), one gets 
as('^r) — 0.334 ± 0.011. We keep conservatively the smallest error, i.e. 

as{ml) = 0.334 ±0.014 — > as(M|) = 0.1204 ± 0.0016 . (6) 

The resulting value is in excellent agreement with the direct measurement of as{Mz) at the 
Z peak, providing a very significant experimental verification of asymptotic freedom. 
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Several recent results motivated the reanalysis of from r decays in ref. [1]. On the 
experimental side, new BABAR measurements of the e+e" annihilation cross section into 
KKn using the radiative return method increase the accuracy on the vector /axial- vector 
fractions in the corresponding r decays, and better results are available on r decays into 
strange final states from BABAR and Belle. On the theory side, the fourth-order term 
in the perturbative expansion of the Adler function was recently calculated [2]. Unitairty 
yields the vector/axial tau hadronic widths in terms of the experimentally accessible spectral 
functions (imaginary parts of two-current correlators 11). The integral can be reexpressed 
using analyticity and Operator Product Expansion (OPE) as the sum of contributions from 
non-perturbative condensates and 5*^°^ the (dominant) perturbative contribution involving 
the Adler function D{s) = —sdU/ds along a circular contour |s| = Sq = M^: 



5(0) 



ds 
s\=so s 



1 



so 



D(s), D(s) 



1 

47r2 



as{s) = as{s)/7T on the contour is obtained by applying the Renormalisation Group Equation 
(RGE) known up to four loops. The phase space involved in the r hadronic width leads to 
an integral for 5*^°^ with a so-called pinched kernel, vanishing for s = Sq. We reexamined the 
convergence properties of the perturbative expansions for the r hadronic widths, and the 
ambiguity between the fixed-order (FOPT) and contour-improved (CIPT) approaches for 
determining along the circle of integration and summing up the series (as the difference 
between the two approaches exceeds the other sources of theoretical uncertainty). Starting 
from the value of in the Euclidean region where OPE is known to be valid, CIPT cor- 
responds to solving the RGE along the circle of integration by small steps, whereas FOPT 
corresponds to a Taylor expansion in powers of ^^(m^) and log(s/so). We compared different 
implementations of FOPT, corresponding to keeping further terms derived from RGE along 
the circle of integration. The main difference between FOPT and CIPT was seen to arise 
from the truncation of higher orders perfomed in FOPT after the contour integral has been 
computed. We identified specific consistency issues for FOPT, which do not exist in CIPT: 
the convergence of the Adler function near the positive axis cut is problematic for FOPT 
(less visible once the Adler function convoluted with pinched weights vanishing at s = m^), 
and the scale dependence of 5^^^ is much more pronounced than for CIPT. 

Possible violations of quark-hadron duality at the r mass scale have been considered using 
specific but crude models (equally-spaced resonances or instanton model). Their effect has 
been found to be well within our quoted overall theoretical uncertainty when we considered 
the V + A spectral function. Due to the very simple models considered, we did not introduce 
additional theoretical uncertainties. We did not investigate the separate vector and axial 
channels, as discussed in refs. [3] and [4], since the more inclusive quantity V + A provided a 
more stable determination of as than V or A separately. This was obtained by performing a 

Rt,v/a and hadronic spectral moments from 



combined fit of the r hadronic width -R^°y/^ 



the ALEPH collaboration of the form R^yjji^ 



!fds (l-4j(;-J^with (fc,0 
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(1, 0), (1, 1), (1, 2), (1, 3), expressed through OPE in terms of 4 unknown quantities D = 4 
gluon condensate, D = 6 and D = 8 effective condensates). The D = 4 quark condensate 
was expressed in terms of quark masses through the Gell-Mann-Oakes-Renner relation, and 
D > 10 contributions neglected, even for higher moments (this approximation may feed back 
in the D = 6,8 condensates). The fit for V+A was indeed more satisfactory than the separate 
V and A fits, yielding ^^(r) = 0.344 ±0.005exp±0.007theo5 consistent with the previous value 
obtained for three known orders, with a 20% reduced theoretical uncertainty. The evolved r 
resuh at the Mz scale was: a^iMz) = 0.1212±0.0005exp±0.0008thco±0.0005evoi, which agrees 
remarkably well with the value directly extracted at this scale (the FOPT value evolved at 
Mz would be lower by 0.0028). 

After the completion of this work, it was pointed out in ref. [4] that the spectral functions 
from the ALEPH data seemed much more correlated than the currently available correlation 
matrix, which indeed does not take into account properly the correlations due to unfolding. A 
study is currently under way to assess these uncertainties precisely, but prehminary studies 
indicate that they have only a limited impact on the correlations between the pinched- 
weight moments used in the analysis of ref. [T]. Another issue was raised in ref. [3], where 
a model was introduced for the higher orders of the Adler function in perturbation theory, 
including the known properties derived from renormalon calculus. This model assumed 
that the perturbative series was saturated by three renormalons (two infrared ones and an 
ultraviolet one) and that the asymptotic behaviour for these renormalons could already be 
used with a good accuracy at order a^. The Borel sum of the model (assumed to yield 
the "true" value of the Adler function) was used to compute S^^^ and compared with CIPT 
and FOPT results at increasing orders of perturbation theory, leading to favouring FOPT. 
In ref. [6], we discussed some assumptions of this model, showing that different conclusions, 
favouring either FOPT, CIPT or neither could be drawn from rather simple variations of this 
model (asymptotic behaviour setting in than at a higher order than a^, presence of several 
subleading infrared renormalons), as well as they depend on the particular weight chosen 
for the analysis. An alternative model [7] led to the conclusion that a different treatment 
of FOPT and CIPT could lead to the same result at high orders for some classes of models 
combining renormalon calculus and conformal mapping. 
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The inclusive r hadronic width provides one of the most precise determinations of the strong 
coupling, since despite the low scale, non-perturbative condensate contributions are below 5% 
of the perturbative correction [1] . Subtracting the condensate terms (PC) from the ALEPH 
data, the experimental value of the perturbative QCD correction relative to the parton model 
decay is (5phen ~ 0.2037 ± 0.0040exp ± 0.0037pc. Theoretically, the perturbative expansion is 
obtained from the vector and axial-vector current correlation functions, whose expansion is 
known to O(a^) [2], by integrating the Adler function along a circle in the complex plane. 
Depending on whether one performs a strict expansion of this integral in ^^(mT-) (fixed order, 
FO) or integrates the running coupling along the circle exactly (contour- improved, CI), one 
obtains 

oo n oo 

n=l k=l n=l 

respectively. In both cases, the dynamical input resides only in the expansion coefficients 
Cn^i of the Adler function, but when one truncates the expansions at some value of n, the 
resulting value of S^^^ is different, since the two expressions treat formally higher-order terms 
differently. The numerical effect is significant, as can be seen from 

4o = 0.1082 + 0.0609 + 0.0334 + 0.0174 + 0.0088 = 0.2288 

4°i = 0.1479 + 0.0297 + 0.0122 + 0.0086 + 0.0038 = 0.2021 (2) 

(using the estimate Cs^i = 283 and ^^(mT-) = 0.34), and increases as more terms are added, 
which points to a systematic problem in one of the two approaches. As will be explained 
below, the fixed-order approach is strongly preferred. Repeating the fixed-order analysis of 
[3] with S^^^^ as given above we obtain 

asinir) = 0.3199 ± 0.0034exp ± 0.0031pc ± 0.0026e5,i tofoll (scale), (3) 

which implies a,(Mz) = 0.1185 ± 0.0004exp 1 aooos (th) ± 0.0002evoi. 

The CI method is often argued to be the method of choice since it exhibits a faster 
apparent convergence and since the expansion of the running coupling on the circle is barely 
convergent near the Minkowskian axis < 0. However, the expansion of the Adler function 
has zero radius of convergence (due to renormalon and other singularities in the Borel plane 
[1]), and is only asymptotic. The finite radius of convergence of the coupling expansion along 
the circle is therefore of subordinate significance - the assumption made in the CI approach 
is that the higher order terms in the Adler function series are negligible in comparison to 
running coupling effects despite their factorial divergence. 
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In order to clarify the discrepancy between the FO and CI approaches, in |3] we included 
in addition to the known c„ i all available information on the structure of the large-order 
behaviour of the Adler function. The analysis is based on the assumption that the Adler 
function series is sufficiently regular, as seen in the known terms, and that sufficiently many 
low-order terms are known, so that one can smoothly connect them with an ansatz for the 
Borel transform that contains the leading singularities. (This does not imply that the series 
is in the asymptotic regime for n ~ 5 - the known terms do not exhibit factorial behaviour. A 
pattern with several singularities being relevant at intermediate n is expected by comparison 
with the exact all-order series in the large- /3o expansion [S].) With such an ansatz for the 
Borel transform we find that only the FO series approaches the "true" result represented 
by Borel sum within its ambiguity, shown as the grey band in the figure above (plotted for 
as{mr) = 0.34), while the CI approximation is systematically too low. In particular, this 
conclusion applies to n = 4, 5, up to which the series is already known (or estimated). 

How general is this conclusion? An important point, specific to the hadronic r width, 
is that the gluon condensate contribution is strongly suppressed. This is refiected in large 
cancellations in the large-order behaviour of the series coefficients Cn,i+gn of the expansion of 
the r hadronic width associated with the first infrared renormalon: {cn,i + gn)/cn,i ~ 
These cancellations are destroyed in the CI expansion, which sums the running coupling 
terms but drops the Adler function coefficients c„.i in higher orders. The CI method 
therefore fails whenever the first infrared renormalon makes the dominant contribution to 
the series of the hadronic r width in intermediate orders. We find that this is always the 
case unless we choose an ansatz in which we suppress this contribution by hand, or allow 
large and unnatural cancellations between the first and second infrared renormalon. 
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One of the largest theoretical uncertainties assigned to the strong coupling constant 
as as determined from hadronic tau decays [2|3|^f5] (for recent reviews see [6l[7f8|9l[T0l[TT] ) 
stems from the differences in the results for Fixed Order Perturbation Theory (FOPT), 
Contour Improved Perturbation Theory (CIPT) and Renormalon Chain Perturbation Theory 
(RCPT). The results of [I] are presented here where it is shown that FOPT and the FOPT 
part of RCPT have much larger errors than usually assumed and agree within their respective 
errors with CIPT. For the non-strange decay width of the r the perturbative part is given 

by 




with the known coefficients Kn=o.A |12|13|14|15|16|17j . The fifth-order term has been es- 
timated to K5 ~ 275 in [17], but the large deviation of the exact K4 from it's prediction 
suggests that a 100% error on is realistic and = 400 ±400 is used in [1]. The methods 
FOPT and CIPT [5] differ in the way ([H) is calculated. In the CIPT approach the /3-function 
is used to numerically solve as{—s) in the complex s-plane by starting with Q!s(m^). The in- 
tegrand is thus calculated in small steps on the circle |s| = and the sum of all pieces gives 
the total integral. For the FOPT method the /3-function and its derivatives are Taylor ex- 
panded in s around Sq = which leads to a power series representation of as{—s) in powers 
of Q;s(m^). The series is truncated at the desired order (here the 5th) in the strong coupling 
and inserted in the integral which becomes solveable now. For the purpose of showing the 
impact of the arbitrary choice to develop around sq = this approach can be generalized 
to first evolve as(m^) to as^ml exp{i(po)) with the numerically solved /9- function and derive 
the Taylor series of (5pert around this new point. In order to preserve the symmetry 

Q;s(m^exp(-z(y9)) = ^^(m^ exp(i(/p))* (2) 

the upper semi-circle is developed around Sq = m^exp(z(y9o) while the lower semi-circle is 
developed around Sq = exp{—iipo) . Other schemes to preserve this symmetry and to 
avoid imaginary contributions to 6pert lead to similar results and are not further discussed. 
Numerically the central result reads [1]: 

5pcrt = a + 0.849 b + (-4.5 Lpo + 0.808) ab+ (1.910 ipo + 5.202) (a^ - b^) + 
(-5.063 Lfl + 5.214 ifio + 26.37) (a^ - 3 a 6^) + (4.297 ifl + 27 Al Lpo + 12.36) (6^ -3a^b) + 
(-9.669 v^o - 101.5 ifil - 71.63 ^0 + 127.1) (a^ -6a^b^ + 6^) + (3) 
(45.56 Lfl - 100.9 ifil - 918.6 ^0 - 521.1) (aH - a b^) + 

(25.63 if^ - 92.90 ip^ - 1221 ipl - 1273 ipo + + 307.8) (a^ - 10 6^ ^ 5 a 6^) + 
(21.76 if^ + 324.8 ifl + 271.8 ifl - 1612 ipo + 0.849 K5 - 1414) {b^ + 5a^b-10a^ b^), 
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with (fio £ [~7r,0], and as{Tn'^exp{i(pQ))/7r = a + ib, which resembles the usual FOPT result 
for (po = and b = 0. 

This FOPT result depends largely on the choice of (fo. FOPT agrees with CIPT around 
(^0 — — 1 but spans over a much larger range of 6pcrt values. Compared to the uncertainty 
from the neglected higher orders this intrinsic error is 4 times larger as none of the choices 
for (fQ should be excluded. The default choice of ipo = leads to the largest possible value 
of (5pert and therefore as from FOPT used to be smaller than from CIPT. The deviation can 
however not be attributed to higher order terms in the series of Sp^rt- 

For the Renormalon Chain Perturbation Theory (RCPT) Spcrt can be written as 

rRCPT _ e _ rFOPT I srFOPT ( a\ 

"pert ~ "renormalon "largo-/3o "pert ) 

where the three terms in the sum refer to the renormalon chain result, the large-/3o re- 
summed result up to the order used in FOPT, and the FOPT result, respectively. For the 
FOPT correction and the fixed order large-/3o correction the same arbitrariness of the choice 
of ipo as discussed above exists, as long as both the FOPT term and the fixed order large-/3o 
term are expanded around the same (po. Therefore the variation of ^pCTt"""" ~~ '^large-^o ^^^^ ¥^0 
is a source of uncertainty in the RCPT approach. The generalized 5ia?gc-^o derived 
from the generalized FOPT solution by setting /3„ = for n > and replacing the Kn with 
/3o" i^n, which are given up to = 4 in [18] and up to n = 12 in [TT]. The numerical 



form is given in [I]. Furthermore the ag in the renormalon part and the fixed order part 
of eq. (jlj) refers to the same quantity. This is probably not the case. The renormalon part 
in [18] is derived from the one-loop coupling in the so-called V scheme, which is 

matched on the 0-loop level to a^^(exp(— 5/3)//^). The problem therefore is that we have a 
coupling constant on the 3-loop level in the FOPT parts, but treat it as a one-loop coupling 
in the renormalon parts. A possible solution would be to use 2-loop matching to go from the 
MS-scheme to the V scheme [19)120] . The generalized RCPT result with 2-loop matching for 

shows a large overlap with the CIPT result and the uncertainties in the deduced strong 
couplings from both theories are much smaller than previously assumed [1]. 

Using the same numerical value for (5port = 0.2042 ±0.0038exp±0.0033non-pert as obtained 
in [To] and used in [11], where the first error is the experimental one, dominated by the non- 
strange hadronic decay ratio of the r, Rry+A and the second is due to the non-perturbative 
and quark- mass corrections, the results for CIPT and generalized FOPT and RCPT read [T] : 

apP'^(m2) = 0.3406 ±0.0047exp±0.0041non-pert± 0.0066x5, 

afOP^K) = 0.3535 ±0.0061exp±0.0053non-pert±0.0208^„ino°?^^, (5) 

af^^^iml) = 0.3440 ± 0.0030exp ± 0.0026non-pert ± 0.0061^„ ± 0.0019^^^., 

All three results agree within the error due to (po which is very large for FOPT but moderate 
in case of RCPT. The difference between CIPT and RCPT is of the same size as the error 
due to (fo for RCPT and the average between both values (and conservatively assigning the 
larger of the two results errors to the average) leads to [1]: 

as{ml) = 0.3423 ± 0.005exp ± O.OOTaa's ± 0.004non-pertlo:oo?^, (6) 

where the fourth error is due to the variation of the renormalization scale. The total theo- 
retical error (including the non-perturbative part) is with ±0.008 only marginally larger 
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than the experimental error. Evolving as given by eq. (j6]) from = 1.7768 GeV to 
mzo = 91.1876 GeV gives [1]: 

a,(m|o) = 0.1213 ± 0.0006exp ± O.OOO8AX5 ± 0.0004non-pcrtlo:ooot ± 0-0002ev, (7) 
where the last error is the evolution uncertainty due to the variation of the thresholds 
niq < mthrcsh < 2mg and the quark masses itself within their respective errors. 
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Current estimates of the size of nonperturbative effects on the precision determination of 
as from hadronic r decays are incomplete for at least three reasons: (i) In none of the existing 
determinations has the effect of duality violations (DVs) been estimated quantitatively. We 
demonstrated in previous work [T], [2] and in this talk that effects from DVs may turn out 
to be numerically significant, and that the use of only doubly-pinched moments of the r 
spectral functions in order to suppress DVs is not reliable, (ii) While the moments used by 
ALEPH[4] and 0PAL[5] in principle probe the OPE up to dimension D=16, the analyses 
retain contributions only up to D=8. This was shown to be not self-consistent [5]. (iii) The 
most precise quoted values for are based on ALEPH data. However, it is now known that 
correlations due to the unfolding of spectral functions were omitted in ALEPH's latest data 
analysis [2]. Our preliminary explorations indicate that the impact of this on the error in as 
is not necessarily negligible. 

The physical reason for the appearance of DVs is the presence of resonances in the vector 
iy) and axial {A) r spectral functions, which neither QCD perturbation theory nor the OPE 
describe adequately. As demonstrated in Ref . [1] , it is possible that the pattern of resonances 
is rather different in vector and axial channels, and this could lead to a significant effect on 
the value of as even if the sum of V and A spectral functions {V + A) looks relatively flat 
in the region between s 1 GeV^ and s = m^. In other words, the assumption that DVs 
can be safely neglected inV + A may turn out not to be justified. Because of this and point 
(ii) above, the estimates of the nonperturbative part of Rr (and other moments) used in the 
literature need to be revisited. 

There is no quantitative theory of DVs in QCD, and in order to investigate this issue in 
more detail, one has to resort to a (physically reasonable) model [1]. This makes it possible 
to explore finite-energy sum rules with simple, nonpinched weights. Figure [T] shows a fit to 
the nonstrange integrated vector spectral function using the sum rule 

t-SQ p j-OO 

I ds pexp{s) = -—. f ds Uope{s) - ds Pdv{s) ■ (1) 

Jo ^^'^ J\s\=so J so 

Here Pexp{s) is the experimental spectral function, available up to sq = m'^, IIope{s) the 
OPE expression for the vacuum polarization, and Pdv{s) the "duality- violating" part of 
the spectral function, i.e., the part not present in Im IIope{s + ie). The contour integral 
over IIope{s) is parametrized by as(m^) and the OPE condensates. The functional form of 
Pdv{s) is not known from first principles. Following Ref. [1], we model it as 

Pdv{s) = K, e'"'" sin {a + (5s) , (2) 

which introduces four new parameters into the fits to experimental data. We restrict sq to 
an interval sq G [smin, '^n^], and we vary Smin in order to check for stability. For the physical 
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Figure 1: Fit to Eq. with the ansatz ^ (red curve) and without the ansatz (green curve). 

motivation of the model of Eq. (j2]), see Ref. [1]. The fit in Fig. [T]was done with OPAL data 
[5], and Smin = 1-5 GeV^. The result for as{m'l) is 

a,(m^) = 0.307(18)(4)(5) (FOPT) , a,(m^) = 0.322(25)(7)(4) (CIPT) , (3) 

where errors are the error from the fit, from varying Smin = 1-5 ± 0.1 GeV^, and from 
varying the estimated coefficient of (as/ir)^ in the Adler function in the range to 566. 

Fig. 1 demonstrates that DVs cannot, in general, be neglected in analyzing hadronic r 
decay data, even integrated versions thereof. This is confirmed by explorations which also 
include the axial channel in the fits. Essentially all existing r-based extractions of as neglect 
DVs, and thus assume model (2) with k, set to zero by hand. Such an assumption is necessar- 
ily accompanied by an uncertainty not accounted for in currently quoted errors for a^. Fig. 
1 shows that it is dangerous to assume, a priori, that this uncertainty can be neglected for 
analyses involving doubly-pinched weights. We are presently studying different combinations 
of V and A spectral-function moments with the goal of quantifying this uncertainty. This 
may help reducing the fitting errors on ag, but at present, if one wishes to be conservative, 
and stick only to results not subject to additional, currently unquantified uncertainties, one 
should treat from r decays as being known to at best the accuracy given in Eq. ([3]). 
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Perturbative Input to Tau Decays 
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The aim of this brief contribution is to recollect the motivation for the evaluation of the 

corrections to the semileptonic tau decay rate, to describe the essential methods and, last 
not least the numerous cross checks which have been performed to establish the correctness 
of the result which, up to now, has not yet been verified by an independent calculation. 

Removing the effect of the CKM matrix element Vud, QED and as well as electroweak 
corrections and subtracting a small nonperturbative piece, the perturbative QCD corrections 
to the Cabbibbo allowed decay rate can be collected in the quantity 6p which has been 
experimentally determined to be close to 0.2 with a relative error of about 2%. In view 
of the relatively large value of as at this low scale, corrections of higher orders play an 
extremely important role in the extraction of the strong coupling constant. In particular it 
seemed important to study the difference between the results based on contour improvement 
and fixed order perturbation theory, that had been observed and studied in [1], and to 
investigate the behavior of this difference upon inclusion of higher orders. 

The key element is the evaluation of the correction to the vector current correlator 
n(g^), more precisely, its absorptive part. Exploiting the renormalization group equation for 
the vacuum polarization function, it can be demonstrated that this requires the evaluation of 
the four-loop contributions to 11, including its constant (non-logarithmic) part and the five- 
loop anomalous dimension of 11, which can be obtained from a suitably chosen combination 
of four- loop propagator integrals [2]. The problem is thus reduced to the evaluation of the 
finite parts of four-loop massless propagators. This is performed in two steps: the algebraic 
reduction to a sum of master integrals f^°'\ multiplied with coefficients C-"^ which depend 
on the amplitude to be evaluated and which are known to be rational functions of the 
space-time dimension d, C = P'^/Q'^, with polynomials and and m + n typically 
of order 60 or larger. (Dimensional regularization is understood throughout.) In [3], [4] a 
method has been proposed which allows to solve for the coefficients C in the limit of large 
d and to represent the coefficients of its (l/d) expansion in terms of Gauss integrals. Given 
sufficiently many terms in the expansion, the rational function can be fully reconstructed, 
and the reduction to masters can be solved purely "mechanically", albeit at the expense 
of enormous algebraic effort. This has required the development of PARFORM [5] and 
TFORM [6], parallel versions of FORM [7] with speedups around ten for twenty processors. 
The analytic results for the 28 master integrals have been obtained [8] using the method of 
"glue and cut" and were obtained in 2004 already. Only recently these master integrals were 
also evaluated numerically [TU], using the program FIESTA [TT], and agreement to four up 
to five significant figures has been observed. 

Using these ingredients, the Adler function and thus the perturbative corrections to 
the tau decay rate and to the R ratio as measured in electron-positron annihilation were 
evaluated in order a^, first for QCD [9] and later for a general gauge theory [14]. A completely 
independent cross check of the reduction to master integrals has not yet been possible. 
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However, recently the four-loop corrections to the Bjorken sum rule have been obtained 
|14] . Using then a generalized version [13] of the Crewther relation [12] which connects in a 
nontrivial way the QCD corrections to the Adler function and to the Bjorken sum rule and 
the QCD beta function, one finds six nontrivial constraints which are indeed fulfilled [13] • 
This observation gives additional confidence in both results. A more detailed discussion of 
various technical aspects can be also found in Ref. [T5] . 

Varying the renormalization scale /i^/M^ in the range 0.4 - 2 one finds [9] af^lM^-) = 
0.322 ±0.004 ±0.02 and af^(M^) = 0.342 ±0.005 ±0.01 for the two options discussed above. 
Here the first error results from the error in (5p, the second one from the fi variation. The 
difference between FO and CI perturbation theory remains unchanged, as anticipated in [1] . 
Taking the mean value between the two results as central value and 0.015 as an estimate of 
the theory uncertainty we find 

a,(M^) = 0.332 ±0.005 ±0.015 (1) 

and, after evolving up to the Z mass 

a,(Mz) = 0.1202 ±0.0019 (2) 

as our final result. This is well consistent with the value of as determined directly from 
the hadronic Z-boson decay rate, including the corresponding NNNLO corrections [T] 

^^(M^) = 0.1190 ±0.0026 (3) 
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Running and Decoupling of ag 

Matthias Steinhauser 
Karlsruhe Institute of Technology 

The /5 function governs the renormalization scale dependence of the strong coupling and 
is known to four- loop accuracy since almost 14 years [1] [2] . There are different ways to solve 
the corresponding differential equation. The prefered method is the numerical solution with 
truncation of I3{as) at the desired order. There are also several approximate (analytical) 
expressions, e.g., the one based on the iterative (perturbative) solution where the result for 
Q!s(/^) is given as an expansion in 1/L = l/ln(/i^/A^) [3]. This formula should be used 
with care, in particular for small renormalization scales yU. If one considers, e.g., /i = Mr 
one observes a shift of +0.004 after including the four-loop corrections and negative shift of 
approximately the same order of magnitude at five-loop level. These numbers have to be 
compared with the current experimental precision which is cited as ±0.005 in Ref. |1] (see 
also the other contributions on as from r decays in these proceedings). 

Next to the running itself also the decoupling of heavy quarks form the running of the 
strong coupling constant is a crucial ingredient of the precision determination of as- Every 
time a heavy quark threshold is crossed one has to apply the decoupling constants which 
relate as with Uf active quark flavours, usually denoted by ai"^'', to the coupling with only 
Uf — 1 active quark flavours. The decoupling constants are obtained by matching rij-flavour 
QCD to the effective theory with the number of quarks equal to — 1. The theoretical 
framework for the calculation of the decoupling constants has been set up in Ref. [5] where 
formulae are given relating /-loop corrections to /-loop vacuum integrals. 



0.1206 
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As a consequence of the decoupling relations as^jj) is not a continuous function of /i 
but has finite steps at the energy scale where the heavy quark is integrated out, //dec- This 
energy is not fixed by theory, should, however, be in the vicinity of the heavy quark mass. On 
general grounds the dependence on /Xdec should become weaker if higher order perturbative 
corrections are included in the analysis. This is demonstrated in the figure above where 
af\Mz) is computed using af\Mr) as a starting point. The decoupling of the charm 
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quark is performed at the fixed scale /^c = 3 GeV and the decouphng scale of the the bottom 
quark fib is varied in the broad range between 1 GeV and 100 GeV. A^-loop running goes 
along with (A^ — l)-loop decoupling. Results are shown for = 1 (upper right dotted line), 
N = 2 (steep dashed line), = 3 (lower dashed line) and = 4 (dash-dotted line). One 
observes a dramatic reduction of the fib dependence with increasing A^ resulting in a quite 
flat four-loop result (Note that the scale on the ordinate only varies by 0.0009.). 

For comparison we show in the figure two more curves which correspond to A^ = 5. They 
incorporate the four-loop decoupling relations [6] [7]. For the unknown five-loop coefficient 
of the P function we have chosen /34 = (solid line) and /S^ = 150 (dashed line parallel to 
the solid one; the normalization corresponding to {/3o, /32, /Ss} ~ {1.92,2.42,2.83,18.85} 
has been chosen). 

From the figure above it is possible to estimate an uncertainty on ai^\Mz) as obtained 
from Q;i^^(M^) due to missing higher order corrections. If we restrict ourselves to a range of 
Hb between 2 GeV and 10 GeV and take the difference between the three- and four-loop curve 
as an estimate for the uncertainty we obtain 6a^s^\Mz) ~ 0.0002. The difference between the 
four- and five-loop (dashed) curve would even lead to 6af\Mz) ^ 0.0003. The variation of 
ai^\Mz) due to the variation of fib leads to an additional uncertainty of dai^\Mz) ~ 0.0002. 
A similar uncertainty is obtained from the variation of fic between 2 GeV and 5 GeV. (This 
can easily be checked with the program RunDec [8].) Thus a total uncertainty of ±0.0004 
(obtained by adding the three uncertainties in quadrature) should be assig ned to JI'\Mz). 
The uncertainties induced by the errors in the quark masses are much smaller. 
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as with GAPP 
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The FORTRAN package GAPP [1] (Global Analysis of Particle Properties) computes 
so-called pseudo-observables and performs least-x^ fits in the MS scheme. Fit parameters 
besides and Mh include the heavy quark masses which are determined from QCD sum rule 
constraints thus affecting and being affected by a^. When possible, analytical expressions 
(or expansions) are used to capture the full dependence on as and the other fit parameters. 

Z-pole observables from LEP 1 and SLC include the Z-width, Tz, hadronic-to-leptonic 
partial Z- width ratios, Ri, and the hadronic peak cross section, (Jhad- These are most sensitive 
to as by far, but the weak angle enters and needs to be known independently. Thus, the 
extracted depends on the set of other, purely electroweak (EW) measurements employed 
in the fits, such as various asymmetries and experiments exploiting parity violation. The 
statistical and systematic experimental correlations of Tz, cxhad and the Re are known, small 
and included. The parametric uncertainties (such as from Mh) are non-Gaussian but treated 
exactly. The theoretical errors in Tz, (Jhad, and the Re are identical, and induce a negligibly 
small uncertainty in Aas{Mz) = ±0.00009, dominated (±0.00007) by the axial- vector singlet 
contribution [2] which is unknown at (9(a^). As in the case of r decays, one may opt for 
either fixed-order perturbation theory (FOPT) or contour-improved perturbation theory 
(CIPT) [3], and we take the difference! as the massless non-singlet uncertainty (±0.00005). 
The ly-width also features a strong dependence, but it is currently not competitive and 
usually interpreted rather as a measurement of a combination of CKM matrix elements. 

The global EW fit excluding r decays (the Z-pole alone) yields as{Mz) = 0.1203 ±0.0027 
(0.1198 ± 0.0028). These results are expected to be stronger affected by physics beyond 
the Standard Model than other as determinations which is the primary reason to include 
another as constraint in the fits as a control. If the new physics affects only the gauge boson 
propagators (oblique corrections) the resulting as{Mz) = 0.1199tQ;oQ3Q hardly changes, while 
allowing new physics corrections to the Z66- vertex gives the lower as{Mz) = 0.1167±0.0038. 

As the aforementioned as control we choose the r lifetime, Tr, not least because of its 
transparent (even if controversial) theory uncertainty. Our master formula |4] reads, 

1 _ KCxpt 

^cxpt ^ ^[ijexpt cxpt 1 _ ^ ^ ^ 291.09 ± 0.48 fs , (1) 

L e,ii 1 directJ ptheo _|_ pthco _|_ pthoo ' \ 

where r^fP^^ = 290.6 (1.0) fs is the directly measured r hfetime [5]. t[BI^p^] = 291.24 (0.55) fs 
is the combination of indirect determinations, using r[Se,^] = ^B'^^^^^ f^^^JT and the experi- 
mental branching ratios, Bf^^^ = 0.1785 (5) and = 0.1736 (5), together with their 13% 
anti-correlation [5]. Decays into net strangeness, S, are plagued by the uncertainty in the 
MS strange mass, rhsirriT-), and a poorly converging QCD series proportional to rn^, so that 
in Eq. ([1]) we employ the measured AS* = —1 branching ratio, Bf^^^ = 0.0286 (7) [5]. 

^This difference has the opposite sign from r decays indicating that their theory errors are uncorrelated. 



37 



The partial r-width into light quarks contains logarithmically enhanced EW corrections, 
S{mr,Mz) = 1.01907 ±0.0003 [6], and reads (employing FOPT as advocated in Ref. [7j), 



^theo 
ud 



Glml\Vud\ 



647r3 
asirrir) 



1 



3 ml 



71 



5M^ 
+ 5.202^ + 26.37^ 



X 



127.1- 



(2) 



a, 

'7l' 



1.393 



5r 



7T 



where Sg collects quark condensate, 5np IH] , as well as heavy and light quark mass effects. The 
dominant experimental and theoretical errors are given in the following tables, respectively: 



source 


uncertainty 


Aas{Mz) 




±0.48 fs 


+0.00039 




±0.0007 


+0.00017 




±0.00022 


+0.00007 


Arrij. 


±0.17 MeV 


+0.00002 


total 




0.00043 



source 


uncertainty 


based on 


Aas{Mz) 


PQCD 


+0.0119 


a^-term 


+0.00167 
-0.00137 


RGE 


/34 = ^579 


m 


+0.00038 
-0.00034 




±0.0038 


m 


+0.00048 


0PE 


±0.0008 


i & [10] 


+0.00012 


total 






+0.00178 
-0.00150 



The perturbative QCD (PQCD) error dominates and is estimated as the af-term in Eq. (jj]). 
It is re-calculated in each call in the fits to access its a^-dependence and features asymmetric. 
It basically covers the range from the higher values favored by CITP down to the lower ones 
one obtains from assuming that the roughly geometric form of FOPT continues. Note that if 
CIPT is used, the error from the renormalization group evolution (RGE) parametrized by the 
unknown 5-loop /9-function coefficient, P^, and part of the PQCD error are correlated. Effects 
breaking the operator product expansion, OPE , are estimated by assuming the instanton 
motivated functional form |9], Aa~^ exp[—27T/as{so)], and adjusting A to the difference 
between the OPE and data curves in Fig. 22 of Ref. [TO]. Our resuh is ajr^] = 0.1174tJ]:g[J}^. 
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Event shapes have been widely used for extracting as values from e"'"e~-collision data. 
They measure the extent to which the energy flow in an event departs from that of a simple 
2-body qq configuration. As as example, the thrust T = equal to 1 for 

a perfect back-to-back two-parton event and 2/3 for a symmetric 3-parton event. Since the 
distribution of thrust values is sensitive to gluon radiation, it can be used to extract a value 
for Os- Various elements go into such an extraction: a perturbative fixed order calculation of 
the distribution (up to NNLO [I]); optionally, a resummation to all orders of logarithmically 
enhanced terms such as a^ln^"~^(l — T) (up to N'^LL for thrust and heavy-jet mass [3], NLL 
for other observables [2]); a model for non-perturbative corrections (either Monte-Carlo based 
or analytical); a choice of the range of event-shape values to be fitted, generally taken to 
ensure that the relation between as and the prediction of the event-shape distribution can 
be considered reliable; and the choice of observables to use. 

A summary of extractions of as is given fig. [H Though they nearly all involve the 
same NNLO perturbative calculation [1], there is a substantial spread in central values. The 
associated uncertainties, nearly always dominated by theory and hadronisation uncertainties, 
are also quite different from one extraction to another. 

For the upper two results of fig. [H the perturbative input is just the NNLO predic- 
tion, supplemented with an estimate of hadronisation corrections derived from the difference 
between parton-level and hadron-level Monte Carlo distributions. The larger theory uncer- 
tainties for the JADE data reflect their lower centre of mass energies, where the larger value 
of as degrades the convergence of the perturbative series. The next two results, which still 
use MC hadronisation, supplement the NNLO calculation with matching to a NLL resum- 
mation, necessary in order to obtain a reliable perturbative prediction near the 2-jet limit. 
The results with resummations tend to be lower than those without (and JADE lower than 
ALEPH) but all remain consistent within total errors. 

The use of hadronisation corrections from Monte Carlo generators brings with it the issue 
that a MC "parton level" is not easily related to the parton level of NNLO or resummed 
calculations. The systematic uncertainty associated with this difference of definitions is 
particularly subtle, insofar as the perturbative calculations implicitly integrate over the non- 
perturbative region. The 3 results that form the central block in fig. [1] attempt to work 
around this issue by using analytical hadronisation models [US], which explicitly attempt 
to remove the double counting between a parametrised hadronisation contribution and the 
non-perturbative component of the NNLO/resummed calculation. Most striking is the SCET 
result, both for its very small error and its apparent inconsistency with the world average. 
In favour of the SCET result, is the finding that as the perturbative order is increased, the 
result remains stable to within estimated uncertainties and the P^r d.o.f. decreases. On 
the other hand, modest variations in the fit range lead to a change in as of ±0.0015, 1.5 
times the quoted theory uncertainty. This hints at missing contributions that are larger 
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Figure 1: Summary of extractions from e"'"e~ event shape and jet-rate measurements. 
Thick (red) error bars are experimental only (statistical and systematic), while the thinner 
(blue) error bars additionally include the theoretical (and hadronisation) systematics. 



than the estimated theory uncertainties. In this context, e.g., it is known that hadron-mass 
effects (and, experimentally, whether or not vr^, etc. are taken to have decayed) can have a 
significant impact on as extractions that use analytical hadronisation models [6] . Studies of 
observables beyond the thrust will also provide valuable insight. 

One possible workaround for the difficulties in treating hadronisation is simply not to 
correct for it and, perhaps, take its whole estimated impact as a systematic error. This was 
essentially the approach for the 5-jet rate extraction, and is potentially an option for jet rates 
because hadronisation tends to be smaller than for event-shape type observables. It may be 
of interest to investigate it also for the 3-jet rate, which seems to have a more convergent 
perturbative series than other observables and has been studied so far with only a subset of 
the e^e~ data. 
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as at NNLO(+NLLA) from (mainly) ALEPH data 
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Event shapes and jet cross sections in e^e~ annihilation were measured to a high accuracy 
at LEP and at previous colliders at lower energies. By comparing the data with the QCD 
description of these observables, it is possible to extract Og- During LEP times, only the next- 
to-leading order (NLO) QCD corrections (plus resummation of logarithmically enhanced 
terms) were available, and the theoretical uncertainty from missing higher order corrections 
was the dominant source of error on the extraction. With the recent calculation of 
next-to-next-to-leading order (NNLO) corrections to three-jet production and related event 
shapes [T], these data can be revisited in view of an improved determination of a^. 

To compare with data, the perturbative NNLO prediction (or matched [2] to next-to- 
leading logarithmic resummation, NNLO+NLLA) is supplemented by quark mass effects (to 
NLO) and has to be corrected for hadronization effects. The standard procedure to estimate 
hadronization effects is to compare parton-level and hadron-level predictions obtained with 
multi-purpose (leading-order, leading-log) event generator programs for the observables. It 
has to be kept in mind that the hadronization models in these generators have been tuned 
extensively to precisely those LEP data sets that are analyzed now for the ag extraction. 

With this procedure, as{Mz) has been determined from ALEPH [3], JADE |5] and 
OPAL [6] data, and the results are summarized in Table [TJ Several important observations 
were made in these studies. First and foremost, inclusion of the NNLO corrections improves 
the mutual consistency of the extractions from the different shape variables. Compared to 
previous NLO studies, the theory error is lowered by almost a factor two. It increases when 
NLLA resummation is included, since this order of resummation is insufficient to match onto 
NNLO. A systematic tension between the ALEPH and OPAL extractions can be observed, 
which can be explained in part by different fit ranges and by the different treatment of quark 
mass effects. 



NNLO 


ALEPH [3] 


0.1240±0.0008(stat)±0.0010(exp)±0.0011(had)±0.0029(th) 


OPAL |6] 


0.1201±0.0008(stat)±0.0013(exp)±0.0010(had)±0.0024(th) 


NNLO 
+NLLA 


ALEPH [4] 


0.1224±0.0009(stat)±0.0009(exp)±0.0012(had)±0.0035(th) 


JADE [5] 


0.1172±0.0006(stat)±0.0040(syst)±0.0030(th) 


OPAL |6] 


0.1189±0.0008(stat)±0.0016(exp)±0.0010(had)±0.0036(th) 



Table 1: Determinations of as{Mz) from event shape data based on NNLO and 
NNLO-I-NLLA, combined with hadronization corrections from multi-purpose event gener- 
ator programs. 



Using hadronization corrections from more recent generators (HERWIG++, combined 
with POWHEG or MC@NLO), results are obtained that are substantially different from 
the PYTHIA values, and outside the band obtained from the different legacy generators. It 
may thus be fair to conclude that the the hadronization models in the generators may be 
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over-tuned in order to compensate the shortcomings of the perturbative description in these 
programs and that the hadronization uncertainty on the as extraction may be larger than 
assumed previously, and could actually be a dominant source of error. In order to circumvent 
this problem, one should either restrict to observables with small hadronization corrections 
(as for example jet rates) or turn to analytic hadronization models which allow a better 
order-by-order matching to the perturbative description of the event shape distributions. In 
two different studies, we pursued both approaches. 

From analytic approaches to hadronization (e.g. dispersive model), it is known that the 
leading power corrections to most event shapes are of order 1/Q. Among the standard set 
of shape variables, this power correction is absent only in (hadronization corrections to 
this variable start at l/Q^). The three-jet rate R^{ycnt)i which derives from Y3 is thus less 
sensitive to hadronization than the event shape distributions. We analyzed the LEPl data 
on i?3 and extracted cng from the measurement at Ucut — 

0.02 as [7] 

as{Mz) = 0.1175 ± 0.0020 (exp) ± 0.0015 (theo). 

Extractions at different jet resolution are highly correlated, and yield consistent results. 

Hadronization corrections to the moments of event shape distributions are additive. The 
moments are therefore very well suited to to disentangle perturbative contributions from 
hadronization effects. After extending the dispersive model for hadronization corrections to 
match onto the perturbative NNLO expression, we used data on event shape moments from 
JADE and OPAL to determine as{Mz) and the hadronization parameter in a combined 
fit [H], yielding: 

as{Mz) = 0.1153 ± 0.0017 (exp) ± 0.0023 (theo). 

The hadronization corrections obtained in this approach are considerably larger than what 
is predicted by the legacy generators. The proper matching of the hadronization corrections 
onto the perturbative expression results in a considerable reduction of the theory uncertainty. 
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Next-to-Next-to-Leading Order and "Classic" Power Corrections 



Bryan Webber 

University of Cambridge, Cavendish Laboratory 
Cambridge CB3 OHE, UK 

Here I report on an analysis of the thrust distribution in e^e~ annihilation, performed 
in collaboration with Richard Davison [1], in which the next-to-leading-logarithmic (NLL) 
resummed prediction of [2] was matched to the next-to-next-to-leading (NNLO) fixed-order 
results of [3], with non-perturbative effects included using the dispersive approach of 

To match the resummed and NNLO fixed-order predictions we use the log-i? matching 
scheme of [2]. The NLL resummed result is expressed as 



R{t) = 



i_t 0-tot dT 



exp[L5(i(asL) + 5-2 («s-^^)] = exp 



00 n+l 
n=l m=n 



where the coefficients G„m are known and L = ln(l/t — 1/tmax + l); ^max = 0.42 being the 
kinematic limit for 5 partons, the maximum at NNLO. The terms with n < 3 are subtracted 
and replaced by the full NNLO result for In R{t). 

To include non-perturbative effects we note that the NLL resummed distribution has the 
form 

1 da g2 

27d 



c 



(2) 



cTtot dT 

where J^{Q^) is the Laplace transform of the jet mass distribution. The idea of |4] is to 
subtract out the infrared contribution to this expression and substitute a non-perturbative 
contribution based on the assumption of a universal low-scale effective strong coupling aes- 
The NNLO perturbative contribution from the infrared region g < /ij is 



S In J^,{Q' 



pert 



71 



In^^ 
Hi 



as (hr) + a1 (hr) — In \- — 



In^ 




(3) 



where hr is the renormalization scale, K2 = C^(67/18 — 7r^/6) — 5nj/9 is the two-loop and 
the three-loop cusp anomalous dimension [5]. At the time of [1] was not known and 

we assumed = which is just equivalent to a change of renormalization scheme. This 

turns out to be a good guess {K3 = 11.0 versus i^l = 11.9 for 5 flavours), although this term 

anyway has a negligible effect. 

Assuming a universal low-scale coupling a^siQ), the non-perturbative contribution from 

the region q < fj,i, where HiuQ -C 1, i.e. 1 — T ^ f^i/Q, is 



S\nJ',{Q' 



n.p. 



2Cf 



71 



Ml 



dqacsiq) i^Q 



2Cf Hi 



(4) 
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The change in the thrust distribution is given by the difference between (jlj) and ([3]), which 
we see is equivalent to a shift by an amount 5T that is a power correction proportional to 
1/Q, modulo some logarithmic dependence if we take hr = Q. Notice that the perturbative 
part ([3]) is the start of a divergent series in asifiR), representing the presence of an infrared 
renormalon. Thus universality of a^s implies that the power correction should decrease as 
we include more terms in the perturbation series. In the fit of [T], the NNLO (order-Og) 
term does represent a significant decrease in the power correction, of 25% relative to NLO. 

The resulting predictions have two free parameters, as(/^ij) and ao(/ii). They were fitted 
to all the available e^e~ thrust data in the cm. energy range 14 < Q < 207 GeV. The 
fitting range was maxj/ij/Q, 0.05} <t< 0.33, giving a total of 430 data points. The best 
fit was at as(91.2 GeV) = 0.1164, ao{2 GeV) = 0.59, with = 466. In view of the large 
contributions from certain data sets, the combined experimental statistical and systematic 
error was quoted in [1] for a Ax^ of 14.6. For A^^ = 1 the experimental error is reduced by 
a factor of 3.4. On the other hand, the theoretical uncertainty was assessed by varying the 
renormalization scale /ir by a factor of \/2 instead of the more usual factor of 2, which would 
double the theoretical error. With these revised error estimates, and taking into account 
loop corrections which give the 'Milan factor' [6], ao — )■ 2J^ao/7i = 0.95 ao, the resulting 
parameter determinations were 

as(91.2 GeV) = 0.11641°:°°^^ , ao(2 GeV) = 0.62 ± 0.02 . (5) 

It is a pleasure to thank the organizers for a very stimulating meeting and to acknowledge 
the support of a Leverhulme Emeritus Fellowship. 
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We perform a global fit of e^e~ data for the thrust event shape at energies Q between 35 
and 207 GeV pQ . We will consider the results with six stages of theoretical improvements: 

1. Matrix elements and fixed order nonsingular terms at order using results from [2]. 
This includes the non-logarithmic term in the hard function at three loops. (A fixed 
order analysis at O(a^) was carried out in [3] using ALEPH data.) 

2. Resummation of the most singular logs with N^LL accuracy using Soft-Collinear Ef- 
fective Theory (SCET) ^. (A purely perturbative ^^hh + O^al) fit to ALEPH and 
OPAL data can be found in [4].) 

3. Profile functions (r-dependent scales /ij, /x^, i?, that account for the multijet 
boundary condition, which ensures fixed-order results are used in the far tail region. 

4. Description of nonperturbative effects with field theory and a fit to a single nonpertur- 
bative matrix element of Wilson lines VLi. ill's uncertainty dominates other OPE terms 
in the tail region, and its Wilson coefficient is included with N^LL + 0(a^) accuracy. 

5. Switching from MS to an R-scheme to define VLi. This ensures Vti and the perturbative 
cross-section are free of (9(Aqcd) renormalon ambiguities. An RGE is used to sum 
large logarithms in the perturbative renormalon subtractions. The fit gives Vti to 16%. 

6. QED final state corrections at 0{q) and NNLL (counting a ~ af); bottom mass 
corrections are included using a factorization theorem with log resummation; 0{a1) 
axial-singlet terms arising from the large top-bottom mass splitting. 

We fit simultaneously for agimz) and VLi in the tail region of the thrust distributions (our 
default 487 bin dataset has r = 1 — T = 6GeV/(5 to 0.33). For fixed Q there is a strong 
degeneracy between as{mz) and Qi. It is lifted by the global fit with different Q values. 
To estimate perturbative uncertainties we performed a 500 point random scan over 12 pa- 
rameters, refitting for each point. This accounts for theory uncertainty from higher orders 
via varying scales, statistical theory errors, and three unknown perturbative coefficients. We 
observe very good convergence when increasing the orders in perturbation theory [1], and 
achieve x^/dof = 0.91 at stage 6. The results are shown in Fig. 1, including a comparison 
with [3] and [1] at the appropriate stages. For the global fit at stage 2 our random scan uses 
the fii variations of [1]. Stage 3 yields a result midway between 1 and 2. 

The most dramatic effect is the —7.5% decrease to ^^(m^) from adding the Qi nonper- 
turbative fit parameter at stage 4. The size of this shift is consistent with a back of the 
envelope estimate based on data [1]. The field theory treatment is in sharp contrast with 
the traditional method of estimating nonperturbative effects employing Monte Carlo gener- 
ators (MC) (used also in [3] and [4]). MC has less perturbative accuracy than level 1, and 
its tune parameters unavoidable include both nonperturbative and higher order perturba- 
tive effects, leading to a double counting in the MC estimate. At the Z-scale the Pythia 
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(Xsimz) from global thrust fits 
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All errors: aAmz) = 0.1 135 ± 0.0012 



perturbative error 
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+ b-mass & QED 
0.1135 ±0.0009 



O.llOL 



Figure 1: Global fit results (solid error bars) at various theoretical stages from [T]. r^i is fit 
at stages > 4. Dashed error bars are ALEPH Q = mz fits from [2] and ALEPH+OPAL 
fits from [4] at stage 1, 2 respectively (hadronization corrections from Monte Carlo are quite 
small and their inclusion therefore does not change the identification with these stages). 



tune has a very small power correction, inconsistent with our determination, and hence 
cannot be used for high precision as determinations. In [5] a global fit was performed at 
NLL respectively + C(af ) with the effective coupling model to treat nonperturbative effects. 
They find as(m^) = 0.1164 ± 0.0022 exp ± 0.0017pert consistent with our stage 2, 4, 5, 6 fits. 
This model agrees with the parametrics of the OPE in the tail region and has renormalon 
subtractions. Advantage of our approach are that the power correction is defined as a field 
theory matrix element, one can include subleading OPE terms, and sum large logs in the 
subtractions. Our final result with all sources of uncertainty is shown in the box in Fig. 1. 

Given the desired precision, an interesting issue is to assess the most appropriate way of 
computing the binned theory distribution. Two approaches being used are: 

, (b) : S(r2,/Xi(r2)) - S(ri,/ii(ri)) , (1) 

where S(r, /i) = (1/cr) /J'dr' dcr/dr(r', /ij(r)) and /ij are jet and soft scales depending on r. 
They are equivalent at a given order in resummed perturbation theory, but differ by notice- 
able amounts from included higher order terms [1]. In particular, S gives an uncertainty 
estimate that is designed for comparison to cumulant data, and (b) is sensitive to contribu- 
tions from r below the bin [ri, 72]. The uncertainties in da/dr are designed for distributions 
and hence (a) is best suited for the bins, (a) was used in [1], while (b) was used in classical 
resummation analyses and [3]. At stage 2 this issue is small for thrust, and only becomes 
noticeable at stage 3. For the Heavy Jet Mass (HJM) event shape it is already noticeable at 
stage 2. The stage 2 thrust and HJM results in Refs. [3] use (b) and yield: 0.1172 vs 0.1220. 
Using (a) yields consistent as{mz) values for thrust and HJM: 0.1169 vs 0.1175, resolving 
the apparent discrepancy in central values. 
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as from five-jet observables at LEP 

Giulia Zanderighi 

Rudolf Peierls Centre for Theoretical Physics, 1 Keble Road, 0X1 3PN, Oxford, UK 

Jet observables represent powerful tools to measure the strength of the strong inter- 
actions, encoded in as- The more jets are present in the final state, the stronger is the 
dependence of the relevant observables on a^. At LEP the highest exclusive multiplicity 
measured is five. However, since five-jet events could be described only at leading order 
(LO) in QCD, the theoretical uncertainty was too large to extract a competitive value for 
as- We recently computed the next-to-leading (NLO) QCD corrections to five-jet events [1] 
using D-dimensional unitarity for the one-loop calculation [2] and Madgraph and MadFKS 
for the real radiation, subtraction, and phase-space integration [3] . We consider two observ- 
ables: the five-jet resolution parameter, 1/45, and the five-jet rate, R^,. Their perturbative 
distributions start at 0{al), and are therefore very sensitive to the value of a^. We find that, 
once NLO corrections are included, ALEPH data is described well by perturbative QCD if 
the world-average of as is used as input of the theoretical prediction. Furthermore, theoreti- 
cal uncertainties are reduced from [—30%, +45%] at LO to [—20%, +25%] at NLO. One can 
then turn this around and use the NLO calculation for a novel extraction of using ALEPH 
data. Before this can be done, one needs to address how non-perturbative hadronization cor- 
rections can be accounted for. Traditionally, they are extracted using parton showers such 
as Herwig of Pythia. We find however that, if this is done, the hadronization corrections 
are very large (~ 50%) and generator dependent. We therefore use the Sherpa event genera- 
tor [3] with its default cluster model [5] to extract hadronization effects. Sherpa implements 
the five-jet LO matrix element exactly, therefore hadronization corrections are less contam- 
inated by missing perturbative effects. When this is done, we find smaller hadronization 
corrections (< 20%). We also observe that while the hadronization corrections per se are 
not negligible, they have a very small effect on the extraction of as at LEP I, therefore we 
neglect hadronization corrections at LEP II. 

In order to extract as we consider each bin of 1/45 and -R5 at a given energy as an observable 
Oi = [Xi,af^^,af'^^] that can be used to measure by solving the equation Tji^j = Ei, 
where Tj is the perturbative theoretical prediction. Hi is the hadronization correction and 
Ei is the experimental measurement. From each bin i we obtain a central value of with 
corresponding errors 

ai = a\ ± (5a*/*"* ± 5a'f^'' ± da'f^^' ± 5af^'^\ (1) 

The systematic and statistical errors are obtained by solving the same equation with Ei = 
Xi ± (j^^^'^^/^y^^ ^ fD^Q perturbative error is obtained by varying the default scale (/i/j = 0.3Mz) 
by a factor two up and down, and the hadronization error is obtained by solving the same 
equation using the Lund hadronization model [B] . We take the fit range as large as possible, 
where our computation is reliable and data are good enough. We then vary the range 
to estimate a fit range error. Full details about the fit ranges and the fitting procedure 
are given in [1]. The result of this procedure it a set of values of as that need to be 
combined. To do this we define a covariance matrix as the sum of individual covariance 
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Figure 1: Breakdown of the errors on as at LEP I and LEP IL 

matrices V = V*^*^* + 1/*^^"* + V'"'^^'' + V^^"^' with 

V,, = 5,, {6al)' + (1 - 6,,) min {5ai)'] , (2) 

where Cij describes the correlation. We assumed that statistical errors are uncorrelated 
at different energies. At a given energy we take 2/45 bins as uncorrelated, R5 bins as fully 
correlated and 1/45 and R5{ycut) to be correlated for ?/cut < 2/45- Systematic errors are taken 
as fully correlated at fixed energy, but uncorrelated at different energies. Perturbative errors 
are taken to be fully correlated for all observables and energies, except for the LEP I/LEP II 
correlation that we neglect. Finally, hadronization errors are assumed to be fully correlated. 
We then compute the weights Wi and use them to obtain the estimate of the average of the 
strong coupling constant and of its error 

N N N N 

Wi = ^(^~^)iiy/ ^ (y~'^)kl, = X] ' "^^("s) = ^ WiVijWj . (3) 

j=l k,l=l i=l i,j=^ 

We obtain the value of as from five-jets (a breakdown of the errors is displayed in Fig. [1]) 

as{M,) = O.U56t'oS- (4) 
This value is compatible with the current world average, but on its lower side. 
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Jet and Event Shape Observables at LHC 

Klaus Rabbertz 
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During the first, very successful year 2010 of LHC operation, which saw an increase 
in instantaneous luminosity by a factor of ~ 10000, an excellent performance of the LHC 
experiments ATLAS and CMS has been observed. The achieved understanding in detector 
operation and behaviour has lead to uncertainties of the order of 11% on the luminosity, 
3-5% on the jet energy scale and about 5-20% on the jet energy resolution depending 
on experiment and phase space. With respect to their impact on QCD measurements three 
general analysis strategies can be differentiated: Absolute measurements like the inclusive jet 
Pt, shape measurements like jet angular variables and event shapes, or ratios of cross sections. 
The latter two exhibit much smaller experimental uncertainties due to cancellation effects 
and are, at least in the beginning, better suited for precise QCD studies in order to constrain 
parton distribution functions (PDFs) and/or determine the strong coupling constant Og- The 
total experimental uncertainties on the jet cross sections currently are of the order of 20 to 
40% for transverse momenta larger than 50 GeV while theoretical uncertainties from the 
PDFs, the renormalization and factorization scales, and non-perturbative corrections amount 
to 8 to 15%. With further improvements and much more data to come the jet cross sections 
are very promising for the future. In the following, some aspects of a shape measurement, 
the azimuthal decorrelation A$jj, and the 3-jet ratio R32 will be discussed. 

Azimuthal decorrelations A$jj, recently published for LHC data by CMS and ATLAS [1], 
are normalized to the total dijet cross section eliminating luminosity uncertainties. They are 
a direct measure of additional activity like a radiated third jet in comparison to a normal dijet 
event perfectly balanced in transverse momentum. Restricting the observation to differences 
in azimuthal angle between the leading two jets reduces significantly uncertainties due to 
the jet energy scale and resolutions. In case of balanced dijet events the azimuthal angles 
between the leading jets are completely correlated to give distances of A$jj = vr. Additional 
radiation of any of the two jets decrease this distance where the extent to which this is 
possible depends on the jet multiplicity as demonstrated e.g. in Fig. 1 of the ATLAS analysis 
in [1]. Comparisons to pQCD at NLO are possible in the range 27r/3 < A$jj < vr with 
NLOJet++ |4] providing the required 3-parton and 4-parton final states. Below 27r/3 mostly 
final states with four or more partons contribute such that the result of NLO Jet++ effectively 
becomes LO only. This can also be seen by an associated increase of the scale uncertainties. 
At A$jj = vr the 3- and 4-parton final states need to be complemented by the 2-parton 
ones requiring 2-loop corrections for a complete NNLO result which is not available as of 
today. This demonstrates that azimuthal decorrelations are a precisely measurable QCD 
observable providing a lot of information on multijet production. Redefining the azimuthal 
distance to A$'- = A$,-,' — tt such that balanced dijets have A$'-- = this quantity can also 
be considered an event shape for which resummed predictions are conceivable [2]. Azimuthal 
decorrelations can also be investigated for leading particles. 

Following a publication of DO, preliminary results by ATLAS and CMS [3] are reported for 
the inclusive 3- to 2-jet ratio R32 versus the scalar sum of all transverse jet momenta Ht- The 
suggested study, however, is far from optimal with respect to comparisons to pQCD. First 
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Figure 1: LO (red) and NLO (blue) predictions including scale uncertainties for an analysis 
similar as in [3] (left) in comparison to the one suggested here (right) requiring a minimal 
Pt of 25% of the leading dijet average < pti,2 >■ 



of all, the LO prediction including the scale uncertainty as shown in Fig. [T]left as upper red 
hatched band is unreliable and gives unphysical values larger than one. The NLO presented 
in blue looks more reasonable and is not far off the CMS data, but the associated scale 
uncertainties are not much smaller with respect to LO. This points to a bad convergence of 
the perturbative series. In fact, ratios -R32 of ~ 0.8 where jets with px > 50 GeV are counted 
even for leading jet momenta oipT,max ~ 1 TeV are not representative for hard processes with 
probabilities related to as at high scales. For that purpose it would be better to require a 
minimal hardness of third jets of e.g. a certain percentage like 25% of the average transverse 
momentum of the two leading jets. An example calculation performed with NLOJet++ in 
the fastNLO framework [3] for such an event selection is presented in Fig. [T] right leading 
to reasonable predictions already at LO and larger reductions of the scale uncertainties at 
NLO at least at high px- The abundant production of jets at highest transverse momenta 
at the LHC will enable many new precision tests of perturbative QCD in the near future. 
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MC tuning with Professor 
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Monte Carlo simulations of high energy physics (and particularly QCD) processes rely 
on phenomenological assumptions for regimes where perturbation theory does not apply. 
Of particular interest for as extractions from event-shape data are hadronisation-effects for 
which two kind of models are available. On the one hand cluster based models, e.g. in the 
codes of Herwig and Sherpa, and the Lund-string model as it is implemented in the Pythia 
event generator. 

Both types have in common that a large number of tunable phenomenological parameters 
are necessary to achieve agreement with data and therefore predictive power. The task of 
tuning such generator parameters can be cumbersome if done by hand due to the high 
dimensionality of the parameter-space and the large number of observables that need to be 
studied for changes due to parameter-shifts. 

The Professor [1] tuning-tool which tries to address the aforementioned problems was 
presented. It is based on the idea of calculating bin-wise parameterisations of the generator 
response to shifts in a certain parameter space and a subsequent numerical minimisation of 
a simple goodness-of-fit measure defined between the parameterisation and corresponding 
experimental data. 

Furthermore, special applications of Professor have been introduced that exploit the 
statistical nature of the tool such as the quantification of tuning uncertainties as well as the 
calculation of so-called "Eigentunes" , a set of parameter points on the /ca-contour of the 
valley that are therefore representative for allowed tuning variations within data- uncertainty 

(Fig. ED. 

Further uses of the parameterisation like calculating the sensitivity of observables to shifts 
in parameter space and an interactive Monte Carlo simulator ( "prof-I" ) have been presented 
as well. 

On request of the organisers, a quick study on a possible a/s dependence of a hadro- 
nisation tune, namely the Pythia6 tune pro-px[T], has been performed. This most recent 
hadronisation tuning was obtained from tuning the generator parameters to event shape 
data on the Z-pole. This is in large parts due to the fact that the development of Professor 
and tunings obtained with it are closely linked to the development of Rivet [2], an application 
designed for generator independent validation by means of built-in analyses that mimic pub- 
lished studies of collider experiments. By the time of the workshop only event shape data 
from the LEP I and II experiments was available in Rivet. An outcome of the workshop was 
the implementation of event shape studies of the PETRA experiments JADE and TASSO 
at energies between 14 and 44 GeVPH] which allowed to study the low-energy behaviour of 
the tune after the workshop. 

Equipped with event shape data and corresponding Rivet analyses at energies between 
14 and 209 GeV pHfS] . the per-bin goodness-of-fit of the tune pYo-p± and data could be 
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Agreement of Pvthia6-tune pro-j'x with data 



Thrust (TASSO) 
Thrust (JADE) 
Thrust (ALEPH) 
■ Yg"'l»" OADE) 

yDurham (ALEPH) 

* C-Parameter (ALEPH) 



100 150 200 

^ [GeV] 

Figure 1: Left: Idea of "EigenTunes" — calculating principal directions of diagonalised 
covariance matrix of a tuning and determination of intersection with ka contour. Middle: 
The 1+ and 1- lines correspond to la Eigentune variations in the positive and negative 
direction of one out of six principal axes. The corresponding parameter points are on the 
la contour of the valley in the vicinity of the central tune point. Right: Goodness- 
of-fit of the Professor hadronisation tune of Pythia6, pro-p^, with data as function of the 
center-of-momentum energy, a/s, indicating a non-universality of the hadronisation tune. 
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studied for the Thrust, C-Parameter and 2-jet-resolution parameter, as a function 

of v^. The result is given in Fig. [H 

The data seems to suggest that Pythia6's hadronisation model or at least the selected 
tuning is not capable of describing the data at all, especially low energies. Further studies 
should include more data, preferably at low energies. It may be fruitful to attempt a retuning 
to all available event shape data simultaneously to check if the available hadronisation models 
are in fact not universal. 

We would like to thank Stefan Kluth for the discussions that lead to the implementation 
of the JADE analysis in Rivet. 
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as from event shapes in e+e : Experimental issues and 

combination of results 

Stefan Kluth 
MPI fiir Physik, Fohringer Ring 6, 80805 Miinchen 
The experimental issues of event shape measurements are summarised in tabled] [T] |2] [3] |1] • 





PETRA/JADE 


LEPl 


LEP2 


Rad. return 


cuts on 4-momentum balance 


n.a. 


cuts on 4-momentum balance, 
cuts on s' 


r pairs 


Nch > 3, reject 1- vs. 3-prong 
events 


Nch>Q 


27 events 


cuts on 4-momentum sum 


4-fermion 
events 


n.a. 


"4-jet event shape" cuts 



Table 1: Summary of experimental effects in e~^e^ event shape analyses. N^h is the multiplic- 
ity of charged particles, 4-fermion events are mainly W^W~ pairs, s' refers to the invariant 
mass of the virtual gauge boson producing the hadronic final state. 



Detector corrections take account of limited acceptance and resolution of the experiments, 
and of QED initial state radiation (ISR) effects, and are performed using MC based unfold- 
ing of the measured distributions. Table [2] gives an overview over the treatment of various 
experimental effects by the LEP experiments, the selection efficients and final background 
fractions, and the resulting experimental systematic uncertainties. The uncertainties are av- 
erage values and can change significantly within a given distribution. The LEPl experimental 
systematics are still to be matched by theory while the LEP2 experimental systematics have 
been reached by recent theory uncertainties. The JADE data have experimental systematics 
varying between 0.7 and 4%. 

Combinations of NLO-t-NLLA as results from the analysis of event shape observables 
were performed by the collaborations using procedures assuming normal distributions for all 
uncertainties. Values of 05 are always evolved to a common scale before combining them. 
The observables are statistically correlated and have correlated systematic uncertainties. 
ALEPH, DELPHI and OPAL use a common method while L3 calculates unweighted aver- 
ages. In the common method experimental uncertainties are taken as partially correlated 
(see below) within an experiment and as uncorrelated between experiments. Hadronisation 
and theory systematics are taken as uncorrelated, and the combination procedures are re- 
peated for different hadronisation models or theory parameters such as the renormalisation 
scale. Statistical correlations within or between observable distributions are measured or 
obtained from MC simulations. For correlations of systematic uncertainties models are em- 
ployed: i) no correlation (p = 0), ii) partial correlation (p = min((Tj, (Tj)^/(crj(jj) or iii) full 
correlation (p = 1). In [6] a combination of all NLO-I-NLLA results from LEP and JADE 
using the common method is shown. However, the recent analyses of LEP or JADE data 
using NNLO(+NLLA) calculations have not yet been combined. 
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ALEPH 


DELPHI 


L3 


OPAL 


Rad. 
return 


find EM 
jets, force 
remaining 
event in 2 
jets 


kin. fit for jet 
and possible 
missing 7 
energies along 
beam 


cuts on 
E ■ N 1 

E^, 

E\ 1 / Eyis 


find isolated 7, force 
remaining event into 4 jets, 
kin. fit for for jet and possible 
missing 7 energies along beam 


4-fermion 
bkg 


force 4 jets, 
cut on 
and 

\mjj — mw\ 


2d cut on Nch 
vs. Bj^ 


force 4 
jets, cuts 

on 7/34, 
■^chi ^clsi 
Ejet 


cut on likelihood from partial 
NLO QCD and EW ME, 2/34, 
Sphericity 


LEPl 
exp. syst. 


< 1% (exp. 
corr.) 


1.6% 


1.5% (exp. 
corr.) 


< 1% (exp. corr.) 


LEPl e 


— 


84.5% 


98.5% ■ A 


88.5% 


LEP2 
exp. syst. 


~ 1% 
(non-rad.) 


~ 1.9% 


1.0-2.5% 
(non-rad 
or 4-f) 


0.8-4% (4-f cuts) 


LEP2 e 




85 - 90% ■ A 


85-90%-A 


~ 80% 


LEP2 fhkg 




5-14% 


20% 


2-6% 


b- Quarks 


NLO 
massive 
ME, no 
corr. syst. 


udsc / udscb 
had. corr., 20% 
syst. 


NLO 

massive 

ME 


neglect, udsc had. corr. syst. 



Table 2: Overview of treatment of experimental effects. E^is is the total measured energy, 
Ncis is the number calorimeter clusters, E^ is the energy of a cluster, E\\ is the energy 
component parallel to the beam direction, au is the angle between jets 3 and 4 after energy 
ordering, Bm is the narrow jet broadening, 2/34 is the jet resolution parameter where the 
event changes from 3 to 4 jets, e is the selection efficiency, fbkg is the background fraction 
and A is the geometric acceptance. 
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Lattice QCD Calculations and ctg 

Andreas S. Kronfeld 
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Lattice gauge theory provides a mathematically rigorous definition of quantum field the- 
ories, such as QCD, encompassing the perturbative and nonperturbative regimes [1]. In 
principle, it is ideally suited for the determination of and the quark masses, which are the 
irreducible parameters of QCD. In practice, the most productive approach to lattice QCD 
has been computational, evaluating the functional integrals of quantum field theory with 
Monte Carlo methods and importance sampling. 

The first such calculations were marred by the omission of sea quarks, often called the 
"quenched approximation." The sea quarks are the most computationally demanding part 
of the calculations and, in the early days, a lack of computer power made the quenched 
approximation a necessary compromise. Quenched QCD calculations are on a footing similar 
to models of nonperturbative QCD. Consequently, the first lattice-QCD determinations of 
required a model-dependent corrections for the omitted sea. In particular, potential models 
could be used to patch up the mass splittings in charmonium and bottomonium [2]. 

Since the start of the new millenium, however, lattice QCD simulations with 2+1 fiavors 
of sea quarks have become routine [3]. The notation "2+1" implies one sea quark tuned to the 
strange quark and two taken with masses as small as possible, for (degenerate) up and down. 
The physical limit is reached by extrapolating numerical data, guided by chiral perturbation 
theory. More recently, the first simulations obviating this extrapolation have been carried 
out [4|5|6] (though not used for — yet). Furthermore, at least two collaborations have 
begun wide-ranging programs of lattice-QCD calculations of 2+1+1 sea quarks — adding the 
charmed sea to the rest pi7] . Thanks to these developments, lattice QCD now plays an 
increasingly important role in many areas, for example, fiavor physics [8]. 

Dating back to Ref. [2], many summaries of perturbative QCD still identify as from lattice 
QCD with the quarkonium splittings. This is outdated. In fact, the possibilities are much 
broader and share many features with determinations of from high-energy scattering or 
the decays of heavy particles. To explain this perspective, let us list the ingredients needed 
to determine a^. One needs 

1. A dimensionless quantity, /?, sensitive to QCD at a (range of) short distance(s), Q'^; 
if R is not dimensionless, one can use Q to make it so; then asymptotic freedom implies 
R = TZ {as{Qs)) + o {{Aqcb / QY ) , where Qs is a specific scale choice for scheme s with 
Qs oc Q, and the power r depends on the observable R. 

2. A theoretical framework — or at least a notion — to separate short-distance scales from 
Aqcd and other long-distance scales. 

3. A perturbative series for the short-distance contribution TZ {as{Qs)), certainly to NLO 
and preferably to higher order, and summing logarithms of scale ratios. 

4. Measurements of R over a range of Q large enough to control the power-law effects, 
denoted (Aqcd/<5)^ in item 1. 
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5. Control of non-QCD physics at scales probed by Q {e.g., electroweak or new physics). 

6. A measurement of Q: usually a calibration, in contrast to event counting for R. 

It is worth emphasizing that the whole procedure can be thought of as dimensionless until 
the last step, where the energy calibration of the calorimeters (high-energy scattering) or 
the r-lepton mass (r decays) is input. 

In lattice QCD, an evaluation of a functional integral replaces the measurement of R. 
One is assured that non-QCD physics does not enter, but, in practice, issues like unphysical 
quark masses and nonzero lattice spacing play an analogous role. We know how to control 
these effects, i.e., how to distill them into an error bar. Many choices of (-R, Q) are possible, 
each with strengths and weaknesses. The dimensionless suite of lattice calculations becomes 
dimensionful when converting from lattice units to GeV — this precisely means deducing the 
scale to which the corresponds (here s labels the scheme chosen in the lattice calculation). 
Since lattice QCD with 2+1 sea quarks reproduces a wide variety of masses, mass splittings, 
and decay constants, the quantity used to convert units is no longer crucial [3], and the error 
on the conversion is sub dominant when propagated to a^. 

The table compiles several recent determinations of from lattice QCD with 2, 2+1, 
or 2+1+1 flavors. Near-term global averages should consider using those set in bold, after 
some dialogue with the authors to ensure a consensus between reviewers and authors of the 
completeness and meaning of the error budgets. The others either are superseded [9IT2] . are 
a re-analysis of a bold entry [TT], omit the possibly large uncertainty from quenching the 
strange sea [TB], or are not yet complete [T7]. Longer-term global averages of will have 
more results from lattice QCD to choose from. 
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Adler function 
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2 Wilson 
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1000 
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m 



* The Adler function's TZ is known to N^LO, but Ref. [14] 's analysis is NNLO. 
This entry gives a^~^\ obtained from Ref. [TBI's tqA^"^'' = 0.62(6) with 
ro = 0.46(1) fm and Eq. (9.5) of the 2006 PDG. 
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HPQCD: as from Lattice QCD 



G. Peter Lepage (for HPQCD) 



Laboratory for Elementary-Particle Physics 
Cornell University, Ithaca, New York 14853, USA 

The HPQCD collaboration has been using lattice QCD simulations for more than a 
decade to determine the QCD coupling from experimental data for hadron masses. That 
history, which has produced consistent results over time, is summarized in the following: 
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HPQCD now uses two very different approaches (red and blue in the plot), with very different 
systematic errors, for extracting the coupling [1]. These give almost identical results for the 
coupling: a^(Mz,n/ = 5) = 0.1184(6) and 0.1183(7). 

The only experimental inputs to a QCD simulation are the small number of well- measured 
hadronic quantities (m^r, rriK ■ ■ ■) used to tune the bare quark masses and bare QCD charge 
in the action. Typically simulations are done for a variety of light-quark masses and lattice 
spacings a (and volumes) in order to extrapolate away lattice artifacts. Extrapolated simula- 
tion results should agree with continuum QCD to within simulation and extrapolation errors 
(of order 1% in many, many tests); there are no free parameters after tuning. We determine 
as by "measuring" various short- distance quantities Y (nonperturbatively) in the simulation, 
and comparing with the corresponding perturbative expansions: Y = Yl'^=i Cn«" (<?!-), where 
the Cn are computed using Feynman perturbation theory. 

Small Wilson loops are the easiest short- distance quantities to measure in a simulation. 
Five-digit precision is routine and the loops are very ultraviolet, and therefore also highly 
perturbative provided (lattice UV cutoff) ir/a is large [2]. HPQCD used 22 different com- 
binations of 8 small loops and the bare coupling, for each of 12 different combinations of 
u-d quark mass and lattice spacing (vr/a = 3.5 — )■ 14 GeV), to extract 22 different values for 
the coupling. Perturbation theory was calculated, through third order, using lattice QCD 
perturbation theory, which includes lattice-spacing artifacts to all orders in the lattice spac- 
ing. Fourth and higher orders were estimated by comparing results from different lattice 
spacings (since q* oc 1/a). Nonperturbative effects, such as condensate contributions, were 
included, and can vary by factors of 100 or more between different quantities. Nevertheless 
all 22 quantities agree to within errors with ^^^^(M^, n/ = 5) =0.1184(6) [I]. 
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To test HPQCD's error estimates, we pretend that perturbation theory is known only to 
first or second order but otherwise use the same analysis, including the estimation of higher- 
order perturbation theory by comparing different lattice spacings. Values for ajjg{Mz) from 
each of the 22 quantities are shown in the following figure: 
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The gray bars show the final result, 0.1184(6), from the full analysis. The first, second 
and third order analyses agree to within their errors of order 0.0092, 0.0018 and 0.0006, 
respectively. Convergence is excellent and the error estimates robust. 

The second approach for determining as used by HPQCD is to measure moments of 
current-current correlators of the form ^^(amoh)^(0|j5(a:, t)j5(0)|0) where j5 = iphl5''Ph and 
/i is a heavy quark (c or b). Low-order moments are short-distance and therefore perturba- 
tive. They also are renormalized quantities, unlike Wilson loops, and so should agree with 
continuum QCD results when extrapolated to zero lattice spacing. HPQCD analyzed four 
moments at approximately 8 different values of the heavy-quark mass between rric and 
and 5 different lattice spacings. The analysis produced very accurate values for the coupling, 
aj^{Mz,nf = 5) = 0.1183(7), and also for the c and b masses [I]. The masses agree with re- 
sults from pure continuum analyses of moments to within about 0.5% — an important check 
on the lattice analysis. 

The agreement between HPQCD's Wilson-loop and current-current correlator determi- 
nations of the coupling is perhaps the most compelling check on each method since the 
systematic errors are very different in the two cases. Many additional tests are discussed 
in p. 
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as from PACS-CS 
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We evaluate the running coupling constant for 2+1 flavors in lattice QCD simulations [1]. 
The Schrodinger Functional (SF) scheme is employed as the intermediate scheme to carry 
out non-perturbative running from the low energy region, where physical scale is introduced 
as 1/Linax; to dccp in the higher energy perturbative region, where a conversion to the MS 
scheme is safely performed: For 3 flavors we have 

Aj| = 2.61192A(^j^, (1) 
^ i,«.)r^exp(-^)exp(-f''.,(J^ + ^-|L)),p) 

where 'g{L) is the running coupling constant in the SF scheme at the box size L, and the 
3-loop beta function f3{g) = —g^ipQ + hig^ + h2g'^) in the SF scheme with 

62 = (0.483(7) - 0.275(5)A^/ + 0.0361(5)A^J - 0.00175(1)A^J) . (3) 

and the universal feo)^i- The SF scheme is the non-perturbative scheme whose scale is 
determined solely by the box size without any other scales. The continuum(a — > 0) limit 
can be taken in this scheme where a is the lattice spacing. 

To convert A^ to o^iMz)-, we first run down at /i = mJmc) using the 4-loop 

beta function and then convert it to at /i = using the 3-loop matching. Next 

we run up a^^{fi) and convert it to a^^{fi) at /i = mb{mb). We finally run up a^^{fi) at 
= Mz{Mz). 

The low energy scale Lmax should be determined by other simulations. We use the scale 
determination in Nf = 2 + 1 QCD simulations at three lattice spacings by CP-PACS/JLQCD 
collaborations [2] whose up-down quark mass covers a rather heavy region corresponding to 
mT,/mp = 0.63 — 0.78. 

Our final results in the continuum limit are 

«,(Mz) = 0.12047(81)(48)(1?73), (4) 

A% = 239(10)(6)(1°,), (5) 

where the first error is statistical, the second one is systematic due to the perturbative 
matching at /i = mc(mc), while the last is due to the continuum extrapolation. 
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Lattice determination of is performed using some perturbative (or short distance) 
quantity calculated non-perturbatively on the lattice. Typical examples are the heavy quark 
potential at short distances and small Wilson loops. One could also consider the vector cur- 
rent correlator in the space-like region. Perturbative expansion for this quantity is available 
to (!) in the continuum theory, and the long-distance contribution is theoretically well 
understood using the Operator Product Expansion (OPE). We propose to use this quantity 
for the determination of [1]. The correlator itself is ultraviolet divergent, but its derivative 
in terms of is finite and unambiguously calculable also in lattice QCD. 

There are potential sources of systematic errors due to (i) discretization effect that could 
be significant at short distances, (ii) truncation error of the perturbative expansion, (iii) 
uncertainty of lattice spacing. These are studied and found to be under good control. Our 
result obtained using our lattices generated with 2-|-l-fiavors of dynamical overlap fermion is 
oi'f'{Mz^ = 0.1181(3)(1;}2) I2]- The error is dominated by the uncertainty of lattice spacing, 
that is significant in our simulation on a relatively small and coarse lattice. 

References 

[1] E. Shintani et al. { JLQCD and TWQCD Collaborations ], Phys. Rev. D79 (2009) 
074510. [arXiv:0807.0556 [hep-lat]]. 

[2] E. Shintani, S. Aoki, H. Fukaya, S. Hashimoto, T. Kaneko, T. Onogi, N. Yamada, 
Phys. Rev. D82 (2010) 074505. [arXiv: 1002.0371 [hep-lat]]. 



64 



as from the ALPHA collaboration 
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The conflict in determining a^jg(yU.): We recall the obvious fact that aMs(/^) °^^y perturba- 
tively defined and not directly measurable. The resulting inherent and practical uncertainty 
vanishes as the renormalization scale becomes very large. In fact the A-parameter is well 
defined from the asymptotic behaviour of aMs(/^ ~^ ^) QCD (in isolation). At low /i, 
where typically the experimental error is small, there is a large theory error in the extraction 
of 

'^Ms(/^) from observables due to the truncation of the divergent perturbative series with 
its intrinsically limited precision. 

Alternatively one would like to connect low and high energies non-perturbatively by 
lattice methods as first suggested for a two-dimensional model pQ. Since on a finite lattice 
that can be simulated, large scale ratios can not be accomodated, low and high energy 
scales are bridged in several steps each changing scale by a factor of 2. This is possible 
by identifying fi = 1/L, where L is the size of the box in which QCD is simulated and by 
using an intermediate (Schrodinger functional) coupling, agp, which is a non-perturbatively 
defined quantity - a finite size effect. 
The master formula of the ALPHA collaboration is 

A^^ 1 r A^i Ai^ 

^ = X ^ X A^^L, X ^ X ^ , (1) 

where A^^/Agp = 2.9065 and typically Lmax ~ 0.5 fm, = 2~'^Ljnax- The first factor 
is computed non-perturbatively and the last three perturbatively. The non-perturbative 
connection of low and high energies yields the relation between asF(l / -^max) and asF(l / Lk) = 
ci;sF(2'^/Lmax) showu in the figure and we then use the perturbative formula giving AgpL^ as 
a function of asF(l/-^fc) in terms of the three-loop /9-function in the Schrodinger functional 
scheme. Since perturbation theory is applied at «sF(l/-^^fc) ~ 0.1 and this renormalization 
scheme has been demonstrated to be well behaved, the perturbative error is negligible. The 
last factor A^/A^ will finally be taken from FT, where we can use the full 4-loop precision 
of [7]. From the numerical results shown by M. Steinhauser at the meeting we estimate 
an uncertainty of 2% on A^/A^ from this step, which is equivalent to decoupling across 

IVL o IVi o 

the b-threshold. In contrast, estimating the ratio A^/A|^ would require matching at the 
charm threshold, which seems very difficult to control perturbatively. We therefore insist on 
the whole strategy being followed with 4 flavours. The approximate decoupling of the charm 
quark at low energy is then implemented non-perturbatively. 

The strategy is transparent and based on minimal assumptions: asymptotic freedom, the 
existence of the continuum limit of the lattice theory and the fact that non-perturbative 
"corrections" vanish quickly as — ?■ oo. Practical issues mainly concern discretization 
effects. The Schrodinger functional is very efficient in rendering simulations with very light 
quarks possible and to achieve the perturbative control (3-loop /3- function) . It does, however, 
introduce unpleasant linear lattice spacing effects. For the desired precision, it is important 
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that the improvement coefficient Ct, which is dominant for removing these hnear effects, is 
known to 2-loops in our discretization and furthermore lattice spacing effects have been seen 
to be small in PT and beyond[2]. 

Status and future: The strategy has been studied for A^f = (pure gauge theory) , and for 
Nf = 2,4 but for the interesting 4-fiavour theory, the ffist factor in eq. ([T]) is not yet know. 



Present results are A^rg 

,(5) 



0.60(5) m, Agro 



0.73(3) (5) 013 compared to A^ro 



0.71(3), Agiro = 0.52(2) from Bethke's world average of 2009, where we use ro 



MS 

0.475 fm. 



For the world average numbers we have not propagated an error in ro = 0.475 fm, since the 
overall error is likely to be revised after this workshop. 

In the near future we will repeat the computation of the scale dependence of asF(/^) on 
faster computers. Smaller lattice spacings and better statistical precision will mean an even 
more precise continuum extrapolation. The future total error of A as we see it now, is given 
in the table and appears to be achievable in a few years. It translates into an absolute error 
of 0.0008 in a,(Mz). 
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Non-pcrturbative running of agp (m) from fi fa IGeV to 
fi ITcV compared to perturbation theory [3]. 
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Error budget. 
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Heavy quarkonium leptonic and non-leptonic inclusive decay rates have historically pro- 
vided ways to extract a^. Ratio of these quantities are very sensitive to ag, if the data are 
sufficiently precise. Determinations of as at the quarkonium masses are important because 
they are among the few ones at a relatively low energy scale. In PDG2006 the reported 
determination of as from T decays was lower than all the other determinations and about 
two standard deviations away the as(M^) average at the time. 

In [1] we have obtained a new extraction of as{Mx(is)) from the decay ratio 

i?^ = r(T(l^)^7X)/r(T(l^)^X) (1) 

X being hadrons. i?^ is calculated in theory using the nonrelativistic QCD (NRQCD) 
factorization for each of the two quarkonium decay amplitudes: 

r = y (^n«MrilS)),mb) (Y(lg)|04fcrmions|^(^g^^ (2) 

where Cn are process-dependent NRQCD Wilson coefficients, rrih is the bottom mass and 
(T(1S')|0^^*'™'°°^|T(1S')) are NRQCD matrix elements of four-fermions operators. They 
contain color singlet contributions that can be calculated in terms of the T(IS') wave function 
and color octet contributions. Adopting the NRQCD counting rules in v {v being the bottom 
quark velocity in the bound state) and as(Mx(i5))/vr ~ ~ a'^{'mf,v), we have included all 
the radiative, relativistic and octet contributions in the expansion of eq.([T]) up to order 
v"^. Higher-order corrections that were not considered are of order v^. We have assumed 
v"^ = 0.08 and rUbV > Aqcd- 

We have used the experimental value R^"^^ = 0.0245 ± 0.0001 ± 0.0013, (where the first 
error is statistical, the second one is systematical) corresponding to the Garcia-Soto (GS) [3] 
parameterization of the r(T(lS') — 7X) data [2], which follows from a QCD calculation. 

Our result is based on: new, precise data from CLEO [2]; a QCD calculation to extrap- 
olate the photon spectrum at low z = [3], based on a combined use of the effective field 
theories NRQCD, pNRQCD and SCET; accurate estimates of the NRQCD color octet ma- 
trix elements in eq.([2]) coming from lattice calculations |1] and from continuum calculations 

By equating the experimental value of to the NRQCD O(f^) theoretical computation 

we get 

as(MT(i5)) = 0.184;H1^ , a,{Mz) = 0.1191°:°°^ , (3) 

very close to the central value of the PDG with competitive errors. Our result is rather 
insensitive to the values of the octet matrix elements because they cancel in the ratio, apart 
from one of them {Os{ ^Si)) which, however, turns out to be small from the lattice calculation 
and does not have a major impact in our results. The determination is valid at next-to- 
leading order in as{mb) and in w^. At this order, terms corresponding to new qualitative 
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features appear (radiative, relativistic, octet corrections), each of them of natural size, but 
whose sum is numerically of order one and hence large. Higher-order corrections 0{v^) are 
expected to be smaller by one order in v since they do not introduce new qualitative features. 
We have used part of the available higher-order corrections to check that this is indeed the 
case [T]. At present, the main uncertainty in our extraction of ag comes from the systematic 
uncertainties in R^^. 

Concerning the extraction of ttg in [5], which due to the small error is dominating the 
final value obtained from heavy quarkonium decays in the PDG2006, there are three main 
differences: (i) On the theoretical side, the color octet NRQCD matrix elements are ignored 
in r(T(lS') — 'jX), whereas we find that they contribute between 30% and 100%. (ii) 
Older data are used, which are fully consistent with, but not as precise as, the more recent 
ones, and an older analysis, which relies on the Field model for extrapolations to low z. 
(iii) The extraction is actually done from r(T(lS') — )■ 7X)/r(T(lS') — l^l^). We believe, 
contrary to a statement in [5], that the latter increases rather than decreases the theoretical 
uncertainties associated to color octet operators. Indeed, whereas the ratio radiative/total 
has the same color octet operators in the numerator and denominator except for one, the 
ratio radiative/leptonic (total/leptonic) has two (three) different color octet operators in the 
numerator and denominator. Furthermore, the leptonic width is known to suffer from large 
higher-order corrections in ag, which introduce further uncertainties. 

Concerning the extraction of the CLEO paper [2], there are two main differences: (i) On 
the theoretical side an old formula was used there, in which the NRQCD color octet operators 
were ignored. This introduces large theoretical uncertainties. In practice, however, we find 
that numerically they are not so important for the final result, (ii) For the total radiative 
width, two numbers are quoted depending on whether the so called Field model or the 
GS parameterization, which is in fact a QCD calculation, are used for the extrapolation 
of the photon spectrum at low z. The final number is given as the average of the two 
procedures. We believe that the use of the Field model, which uses a parton shower Monte 
Carlo technique to incorporate the effects of gluon radiation by the outgoing gluons in the 
decay, introduces an unnecessary model dependence that moves the actual central value and 
artificially increases the errors. Our final results are similar to the ones presented in [2] for 
the GS parameterization. 
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Beyond its intrinsic importance for phenomenology and technological developments, 
knowing the precise value of can provide important information on new physics beyond 
the standard model. Here we discuss two such examples where as may play a crucial role in 
testing new physics scenarios: (1) Grand Unification and (2) the Multiverse. 

Grand Unification in Higher Dimensions — Grand unification provides a beautiful 
understanding of the gauge quantum numbers of the quarks and leptons. It also leads to a 
prediction of as through unification of the three gauge couplings at a high energy, allowing 
us to test the idea. In fact, in supersymmetric grand unified theories, the predicted value 
of as agrees with the experimental value at the 10% level. However, because of unknown 
threshold corrections at the unification scale, the prediction is somewhat model dependent; 
and in typical models, it is around as — 0.130, so that the agreement is not perfect [1]. 
Moreover, theories often suffer from severe phenomenological problems, such as doublet- 
triplet splitting, too fast dimension five proton decay, and unwanted mass relations between 
light generation fermions. 

Grand unification in higher dimensions elegantly solves these problems [2]. In this frame- 
work, the unified gauge symmetry is broken due to a compact extra dimension(s) of order the 
unification scale. The gauge and Higgs fields propagate in the bulk of higher dimensional 
spacetime, leading to automatic mass splitting between the standard model and unified 
gauge bosons, as well as between doublet and triplet Higgs fields. The structure of higher 
dimensional mass terms automatically suppresses dangerous dimension five proton decay. 
Matter fields may either propagate in the bulk or localized on one of the branes, i.e. the 
"edge" of compact space, which can reproduce a complicated pattern of quark and lepton 
masses and mixing without contradicting unification at high energies. A schematic picture 
of this framework is depicted in Fig. [1] 

An interesting aspect of the above framework is that threshold corrections at the unifi- 
cation scale are fixed in terms of only a few parameters. In the minimal model of super- 
symmetric SU{5) in 5 dimensions, it is only the size of the extra dimension M*L (in units 
of the cutoff /string scale) that we need to know to compute the threshold effects. Together 
with the assumption of strong coupling at the cutoff scale, which fixes M*L ^ O(IOO), we 
can predict as without any free parameter. The value obtained is 



where the first and second errors represent uncertainties from the superpartner spectrum and 
M^L, respectively. The former can be reduced/eliminated if superpartners are discovered 
at the LHC. It is interesting that the minimal model indeed improves the prediction over 
conventional 4 dimensional supersymmetric grand unification. 




(1) 
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in the "bulls" 



Figure 1: The basic framework of higher dimensional grand unification. 



Precision Higgs Mass Prediction from the Multiverse — The origin of the elec- 
troweak scale is one of the major mysteries in the standard model. While its smallness 
(compared with the Planck scale) may be ensured by new physics at the TeV scale, it is 
possible that it arises simply as a result of environmental selection in the multiverse (as does 
the cosmological constant). The low energy theory may then be simply the standard model. 
Is there any experimental handle on such a scenario? 

While supersymmetry is no longer needed at the TeV scale in this scenario, it may still 
exist at very high energies, e.g., as a consequence of string theory. Suppose supersymmetry 
exists at some high energy scale rh, i.e. superpartners have masses of order rh. In this 
case, the theory below rh is just the standard model, but with the Higgs quartic coupling 
determined by the weak gauge couplings: \ = {g"^ + 5''^)/8 in the limit where the low energy 
Higgs arises from one Higgs supermultiplet (large tan /3). This allows for predicting the Higgs 
boson mass through renormalization group equation. The result is surprisingly precise: 

V 1-3 GeV J V 0.002 

+ 0.14 GeV (logioY^T^^) + 0-10 GeV ± 0.5 GeV, (2) 

where 5 represents possible threshold corrections, typically of 0(0.01). The uncertainty is 
dominated by experimental errors on the top quark mass and as, which can be reduced. 
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Rehable quantitative predictions for cross sections at the LHC require control over uncer- 
tainties arising from QCD. These errors take two primary forms: our inability to calculate 
to high enough order in the QCD perturbative expansion, and the imprecise determination 
of the non-perturbative input parameters of QCD such as the parton distribution functions 
(PDFs) and the strong coupling constant. Impressive progress has been made recently in cal- 
culating the perturbative cross sections needed in LHC analyses. Fixed-order cross sections 
for high-multiplicity scattering processes are becoming available through next-to-leading or- 
der (NLO) using powerful new calculational techniques [1]. Next-to-next-to-leading order 
(NNLO) calculations of benchmark processes such as Higgs and electroweak gauge boson 
production are available in the form of flexible numerical programs which are easy to incor- 
porate into experimental analyses [2]. For several important LHC processes, the dominant 
uncertainties now come from determinations of PDFs and as. 

The most striking such process is Higgs boson production through gluon fusion. This 
process is famously sensitive to large perturbative QCD corrections. The theoretical commu- 
nity has invested several decades into understanding the gluon-fusion mechanism. A simple 
effective-theory framework for this process exists, and the NNLO QCD corrections, NLO 
electroweak corrections, and the resummation of large threshold logarithms are available [3]. 
Recent predictions for the inclusive cross section used by the LHC collaborations have been 
given in Ref. from which the following Table of -y/i = 7 TeV results has been taken. 



Mh (GeV) 


a(pb) 


Scale (%) 


PDF4LHC (%) 


120 


16.63 


+7.2 - 7.9 


+7.6 - 7.0 


150 


10.52 


+6.6 - 7.4 


+7.6 - 7.5 


180 


6.76 


+6.2 - 7.0 


+7.5 - 7.8 


210 


4.74 


+6.0-6.7 


+7.5 - 7.9 


240 


3.59 


+5.9-6.4 


+7.7-8.0 


300 


2.42 


+5.8-6.0 


+8.0-8.3 


400 


2.03 


+5.9-5.4 


+8.8-8.6 


600 


0.336 


+6.1 - 5.2 


+ 10.1 - 9.4 



The uncertainties arising from PDFs and the strong coupling constant, denoted as 'PDF4LHC' 
in the Table (named for the procedure by which they are obtained, described in Ref. |4]) 
dominate for all Higgs masses over the estimated errors arising from unknown higher-order 
QCD corrections denoted as 'scale.' The situation is potentially worse than indicated, as re- 
cently emphasized [5]. As the Higgs cross section is proportional to a ~ and receives 
large 0{a\) corrections, it is sensitive to the myriad assumptions that enter extractions 
of as and the PDFs. The current PDF4LHC error estimate neglects PDF extractions for 
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which the Higgs production cross section is up to 25% lower than the sets utihzed in the 
Table above. These sets leading to a lower production cross section typically also obtain 
lower values of as{Mz) from their fit of the available data. For example, the ABKM set 
which gives the lowest predictions has as{Mz) = 0.1135 ± 0.0014 [H], while the HERA- 
PDF set which gives the second lowest result has as{Mz) = 0.1145. Both as values ob- 
tained in these PDF extractions are many standard deviations away from the current world 
averagea5(M^) = 0.1184 ± 0.0007 [7]. A proposal to add an additional 'theoretical' uncer- 
tainty on the extraction oias leads to an increased error that partially resolves the differences 
between extractions at low-to-mid Mh [5]. The differences at high Mh are primarily due to 
differing gluons in the various sets, and are not accounted for by this procedure. Neverthe- 
less, further exploration of this issue by the community, and the consideration of whether a 
'theoretical' error on as should be added to the Higgs error budget at the LHC, would be 
welcome. 

One may also ask if may be constrained by measurements in other channels at the 
LHC. Event shapes in jet production may eventually permit an extraction with a 3% uncer- 
tainty, as discussed by Klaus Rabbertz in these proceedings. Studies of -|- jet production 
also appear promising. A conservative estimate of the theoretical error on the ratio of ly -|- 1 
or more jets over the inclusive W production cross section is already at the few-percent 
level [8], with the prospect for further improvement once the NNLO corrections to the + 1 
jet process are known. Luminosity and other systematic errors cancel in the ratio; the statis- 
tical and remaining systematic errors will improve as more data is collected. Further study of 
whether the strong coupling constant can be measured using these observables is warranted. 
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The International Linear CoUider (ILC) [1] is a proposal for an electron positron col- 
lider operating at a centre of mass energy oi Q = 500 GeV (with an upgrade option to 
1 TeV) and an luminosity of 2 • 10^^/cms. We illustrate the precision gain of an mea- 
surement at 500 GeV by the reanalysis of JADE [2] and OPAL [3] event shape measure- 
ments, where as is determined to next-to- next-to-leading order, with the results ^^(m^o) = 
0.1210 ± 0.0061 (^^(m^o) = 0.1172 ± 0.0051) in the JADE analysis without (with) inclu- 
sion of resummed logarithms in the NLLA approximation, and as{Tnzo) = 0.1201 ± 0.0030 
{as{mzo) = 0.1189 ± 0.0041) in the corresponding OPAL analyses. With a precision of 
2.6%-5.0%, these are among the best measurements. In comparison to corresponding NLO 
analyses, the renormalisation scale uncertainty and the scatter from different variables are 
reduced. 

We estimate the uncertainties of a corresponding as{mzo) measurement at 500 GeV as 
follows. The cited luminosity is about an factor 1000 higher than at LEP 1, the hadronic 
cross section is lower by this factor [5]. The selection efficiency is slightly lower than at 
LEP 2 [6]. Therefore a statistical precision of 0.0001 can be expected after taking data for 
a few years. The detector uncertainties in the JADE and OPAL analysis were of the order 
of 0.001. The ILC detector will be very hermetic and supply good tracking and calorimetry, 
and so this order of magnitude is expected to still hold [7j. The uncertainties induced by 
residual background after the JADE or OPAL selection cuts are of the order of 0.001. The 
effects from W and Z pairs rise with centre-of-mass energy. Above ~350 GeV there will 
be additional background from electroweakly decaying top-pairs. The W contribution can 
be suppressed [8] by beam polarisation, and the properties of W and Z pairs have been 
measured properly at LEP. Therefore an uncertainty can be expected which is still of the 
same order [7]. The uncertainties of as{Q) related to the hadronisation correction are of 
the order of 10~^...10~'^ for JADE and OPAL. They fall steeply with centre-of-mass energy 
as can be understood from analytical power correction models. Their expected order in an 
measurement at 500 GeV is 0.0001 after evolving the as value back to 91 GeV. Estimates of 
the scale uncertainties from an NLO and NNLO measurement are given in Tab.[Tl Measuring 

a, (91 GeV) a, (500 GeV) a, (91 GeV), evolved 
measurement estimate from 0^(500 GeV) 



NLO, missing: 
NNLO, missing: a 



4 



0.1192 ±0.0047 0.0959 ±0.0024 0.1192 ± 0.0038 
0.1205 ±0.0027 0.0967 ±0.0011 0.1205 ± 0.0017 



Table 1: as measurements at 91 GeV with the scale uncertainty [3] [4], their evolution to 
500 GeV and estimates of the scale uncertainty at 500 GeV. 

as at 500 GeV instead of 91 GeV reduces the scale uncertainty to 80% in the NLO case, and 
to 60% in the NNLO case. 



73 



To study the measurement of thrust by the SiD detector at the ILC, we process u, d, 
s, c, b - quark pairs, generated by whizard at 500 GeV with PYTHIA parton shower and 
hadronization. The effects of initial state radiation and beamstrahlung are included and a 
cut of m > 475 GeV is imposed on the reconstructed quark pair invariant mass. Comparing 
the thrust distribution from truth particles (with mean (T) = 0.953) and reconstruction 
((T) = 0.952), no systematic deformations are seen. Studying the reconstructed thrust 
values versus the MC values. Fig. [H the migrations in the range T > 0.7, which is measured 
by OPAL [1], are small. 



reconstructed 
1.0-- 

0.9-- 

0.8 

0.7-- 

0.6-- 



0.5 
0.50 



■I 



.: ■ 'lllli 



—\ h- 

0.60 0.65 



+ 



—\ h- 

0.80 0.85 



—\ h- 

0.90 0.95 



0.55 



0.70 



0.75 
MC 



1.00 



Figure 1: Thrust values of the reconstructed particles versus values of the MC input particles. 
We would like to thank Norman Graf for help with the SiD Monte Carlo samples. 
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The Principle of Maximum Conformality 

Stanley J. Brodsky and Leonardo Di Giustino 
SLAG National Accelerator Laboratory, Stanford University 

A key problem in making precise perturbative QGD predictions is the uncertainty in de- 
termining the renormalization scale of the running coupling ^^(yU,^). It is common practice to 
guess a physical scale fi = Q which is of order of a typical momentum transfer Q in the pro- 
cess, and then vary the scale over a range Q/2 and 2Q. This procedure is clearly problematic 
since the resulting fixed-order pQGD prediction will depend on the renormalization scheme, 
and it can even predict negative QGD cross sections at next-to- leading-order [1]. Other 
heuristic methods to set the renormalization scale, such as the "principle of minimal sensi- 
tivity" [2], give unphysical results [3] for jet physics, sum physics into the running coupling 
not associated with renormalization, and violate the transitivity property of the renormaliza- 
tion group [1]. Such scale-setting methods also give incorrect results when applied to Abelian 
QED. Note that the factorization scale in QGD is introduced to match nonperturbative and 
perturbative aspects of the parton distributions in hadrons; it is present even in conformal 
theory and thus is a completely separate issue from renormalization scale setting. 

Scales in QED: There is no ambiguity in setting the renormalization scale in quantum 
electrodynamics: In the standard Gell-Mann-Low scheme for QED, the renormalization scale 
is simply the virtuality of the virtual photon. For example, in electron-muon elastic scatter- 
ing, the renormalization scale is the momentum transfer t; i.e., a{t) = a(to)/(l — n(t, tg)) 
where n(t, tg) = (n(t) — n(to))/(l — n(to)) sums all vacuum contributions in the dressed 
photon propagator, proper and improper. Although the initial choice of renormalization 
scale to is arbitrary, the final scale t is not. In the case of muonic atoms, the modified 
muon-nucleus Goulomb potential is precisely (y{q'^)/q'^. One can use other renormalization 
schemes in QED, such as MS scheme, but the physical result will be the same after allowing 
for the displacement of scales. For example, if >> mj, 0!jjg{e~^^^t) = acM-iit)- The 
same underlying principle for scale setting must hold in QGD since the up terms in the QGD 
(3 function have the same role as the lepton A'^ vacuum polarization contributions in QED. 

PMC and BLM: The purpose of the running coupling in gauge theory is to sum 
all terms involving the /3 function; when the renormalization scale /i is set properly, all 
nonconformal /3 7^ terms in a perturbative expansion arising from renormalization are 
summed into the running coupling. The remaining terms in the perturbative series are 
then identical to that of a conformal theory; i.e., the theory with /3 = 0. The divergent 
"renormalon" series of order a^P^n). does not appear in the conformal series. Thus as in 
QED, the renormalization scale /i is determined unambiguously by the "Principle of Maximal 
Gonformality (PMG)". This is the principle underlying BLM scale setting [5] An important 
feature of PMG is that its QGD predictions are independent of the choice of renormalization 
scheme. The PMG procedure also agrees with QED in the Nc limit. In the case of 
e"'"e~ annihilation to three jets, the BLM/PMG scale is set by the gluon jet virtuality. 

Global PMC Scale: Ideally, as in the BLM method, one should allow for separate 
scales for each skeleton graph; e.g., for to electron-electron scattering, one takes a{t) and 
a{u) for the t-channel and w-channel amplitudes, respectively. Setting separate scales can 
be a challenging task for complicated processes in QGD where there are many final-state 
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particles and thus many possible Lorentz scalars qf. However, one can obtain a useful first 
approximation to the full BLM-PMC scale-setting procedure using a single global scale fi 
which appropriately weights the individual BLM scales. The global scale [6] can be deter- 
mined by varying the subprocess amplitude with respect to each invariant, thus determining 
the coefficients fi oilogq]^ / nl in the amplitude; the global PMC scale is then fi^ = CIli{qf)'^% 
where the weight Wi = fi/ fj. C is the scheme displacement; e.g., C = e~^^^ for MS. 

Commensurate Scale Relations (CSR) |I7|: Relations between observables must be 
independent of the choice of scale and renormalization scheme. CSRs are thus fundamental 
tests of theory, devoid of theoretical conventions. For example, the PMC relates the effective 
charge ag^(Q^) determined by measurements of the Bjorken sum rule, to the effective charge 
aji{s) measured in the total e+e~ annihilation cross section: [1— ag^(Q^)/7r] x [l-|-aij(s*)/7r] = 
1. Because all /3 7^ nonconformal terms are absorbed into the running couplings using PMC, 
one recovers the conformal prediction [8]; in this case, it is the Crewther relation [9]. The 
ratio of PMC scales \fs* /Q — 0.52 is set by physics; it guarantees that each observable 
goes through each quark flavor threshold simultaneously as and s are raised . Thus by 
applying PMC, the conformal commensurate scale relations between observables, such as 
the Crewther relation, become valid for non-conformal QCD at leading twist. 

Conclusions: The PMC provides a consistent method for determining the renormal- 
ization scale in pQCD. The PMC scale-fixed prediction is independent of the choice of 
renormalization scheme, a key requirement of renormalization group invariance. The results 
avoid renormalon resummation and agree with QED scale-setting in the Abelian limit. The 
PMC global scale can be derived efficiently at NLO from basic properties of the PQCD cross 
section. The elimination of the renormalization scheme ambiguity using the PMC will not 
only increases the precision of QCD tests, but it will also increase the sensitivity of colliders 
to new physics beyond the Standard Model. 
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Workshop on ag 2011: Summary 



Siegfried Bethke 
Max-Planck-Institute of Physics, Munich, Germany 

At this workshop, a large number of recent and newest developments in the field of 
precision measurements of have been presented, and many in-depth discussions about 
specific topics which arose before and during this meeting were performed. At the end of 
this workshop, the status can be briefly summarised as follows: 

• A selection of the newest significant results, not included in the 2009 world sum- 
mary (c.f. the Introduction to this workshop, these proceedings), with potential to be 
included in a new world average 2011, gives: 

1. R. Abbate et al., "Thrust at N^LL with power Corrections and a Precision Global 
Fit for a,(Mzo)", arXiv: 1006.3080, Phys. Rev. D83, 074021 (2011): a,(Mzo) = 
0.1135 ± 0.0011. Hadronization power corrections treated with field theoretic 
matrix elements fit to experimental data. 

2. G. Disssertori et al., "Precise determination of the strong coupling constant at 
NNLO in QCD from the 3-jet rate in e+e" annihilation at LEP", arXiv:0910.4283, 
Phys.Rev.Lett 104(2010) 072002: a,(Mzo) = 0.1175 ± 0.0025. 

3. DO Collaboration, V.M. Abazov et al., "Determination of the strong coupling 
constant from the inclusive jet cross section in pp collisions at ^/s = 1.96 GeV", 
arXiv:0911.2710, Phys.Rev.D80 (2009) 111107: a.(Mzo) = 0.1161 ±0.0045 (QCD 
in NLO). 

4. OPAL collaboration, G. Abbienti et al., "Determination of using OPAL hadronic 
event shapes at ^/s = 91 — 209 GeV and resummed NNLO calculations" , subm. 
to Eur.Phys.Jour.C: ^^(Mzo) = 0.1189 ± 00041. 

5. A. Pich, "a^ Determination from r Decays: Theoretical Status" , arXiv: 1001. 0389, 
Acta Phys. Polon. Supp. 3, 165-170 (2010): a,(Mzo) = 0.1213 ± 0.0014. 

• The most imminent problems remaining after all discussion are: 

(a) The new method and result from e~^e~ Thrust distribution (result no. 1 above), 
using full NNLO perturbative QCD plus N^LL power corrections from a field 
theory ansatz provides a rather small value of as(Mzo) and - moreover - a very 
small overall error. This, if included in the 2009 world average (replacing the 
standard e~^e~ results from event shapes; c.f. the "Welcome and Introduction" 
article to this workshop), results in an average of as(Mzo) = 0.1174ib0.0006, with 
a rather bad overall of 17 for 6 d.o.f. This single measurement then deviates 
by more than 4 standard deviation from the exclusive mean (of 0.1185 ± 0.0007), 
indicating that either this determination suffers from a large and so far unknown 
systematic shift, and/or its assigned error is largely underestimated. 

(b) The situation in the sector of ag from r-decays is still unsolved. Even after inten- 
sive discussion throughout this workshop, no agreement on a possible preference 
for one of the different perturbative approaches (FOPT vs. CIPT) was achieved. 
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There was some scepticism on the various results from lattice gauge theory pre- 
sented at this workshop, especially on the reliability of the very small quoted 
uncertainties. It was suggested to include more than the one result which was 
part of the 2009 summary. 

The reliability of the "standard" as determinations from hadronic event shapes 
and jet rates in e^e~ annihilation final states, based on hadronisation corrections 
derived from Monte-Carlo models, was also discussed and criticized. 

General concern was raised about continuously shrinking systematic (and theo- 
retical) uncertainties assigned to some of the results. Defining theoretical uncer- 
tainties by varying the QCD renormalisation scale within a factor of two, up and 
down, of the canonical scale - often the cm. energy of the collision, may be a 
gross under-determination in some cases, especially in case of QCD in NLO only, 
but also in NNLO, where a small dependence on the scale may be accidental, and 
modifications of the renormalisation scheme would also be mandatory. 

A new attempt to extract an updated world average value of as{Mio) will therefore 
have to rely on a careful selection and treatment of the results to be included. As a very 
preliminary first look into the new situation in 2011, the results listed in table 1 are selected 
to calculate a new world average. Compared to the previous study in 2009, results no. 3 and 
4 as listed above have been included, while results no. 1, 2 and 5 were not (yet) accounted 
for due to the reasons given above. The weighted average is as{Mio) = 0.1183 ± 0.0006 with 

= 5.0/8 d.o.f. Assuming an overall correlation factor of 0.37 to adjust to 1/d.o.f. then 
gives as very preliminary world average 2011: 

a^(Mzo) = 0.1183 ±0.0010 . 



Process 


Q [GeV] 


as(Mzu) 


excl. mean as(Mzo) 


std. dev. 


r-decays 


1.78 


0.1197 ±0.0016 


0.11809 ±0.00109 


0.8 


DIS [F2] 


2 - 170 


0.1142 ±0.0023 


0.11866 ±0.00132 


1.7 


DIS [e-p jets] 


6 - 100 


0.1198 ±0.0032 


0.11827 ±0.00097 


0.5 


Lattice QCD 


7.5 


0.1183 ±0.0008 


0.11838 ±0.00164 


0.0 


T decays 


9.46 




0.11832 ± 0.00094 


0.1 


e"'"e~ [jets & shps] 


14 - 44 


0.1172 ±0.0051 


0.11835 ± 0.00094 


0.2 


pp incl. jets 


50 - 145 


0.1161 ±0.0045 


0.11831 ± 0.00097 


0.5 


e"'"e~ [ew prec. data] 


91.2 


0.1193 ±0.0028 


0.11829 ± 0.00095 


0.3 


e"'"e~ [jets & shps] 


91 - 208 


0.1208 ±0.0038 


0.11826 ± 0.00099 


0.7 



Table 1: Preliminary summary of recent measurements of as(M2o) (Feb. 2011). The rightmost two 
columns give the exclusive mean value of as{Mio) calculated without that particular measurement, 
and the number of standard deviations between this measurement and its exclusive mean. 
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